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ABSTRACT 


Ап experimental apparatus has been constructed and 
tested to examine the convective heat transfer in straight 
and curved ducts of rectangular cross-section. The channel 
can be configured with both walls independently heated at a 
constant heat flux, or one wall at a constant heat flux and 
the opposite wall adiabatic. Local and average Nusselt 
numbers can be calculated, and used to evaluate the effect 
of curvature on heat transfer. 

Experiments were conducted at steady state for turbulent 
flow with one wall heated at a constant heat flux and the 
opposite wall adiabatic. The heat transfer characteristics 
of the straight and curved sections, on both the inner and 
outer walls, were compared. The heat transfer rate of the 
concave curved wall proved to be enhanced over both the 


convex curved wall and both straight sections. 
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I. INTRODUCTION 


А. TAYLOR-GORTLER VORTICES 

Since the early part of the twentieth century, 
considerable research has shown that the fully developed 
laminar flow along a concave wall does not remain two- 
dimensional [Refs. 1, 2, 3]. The flow instead forms a 
system of spiral vortices, of counter rotating pairs, whose 
axes are aligned in the principle flow direction. This 
phenomenon is the result of the variations in the 
centrifugal forces acting on the fluid particles, and is 
known as Taylor-Gortler vortices. Figure 1 illustrates the 
type of fluid motion just described. 

In a channel that is curved in the streamwise direction, 
those fluid particles located in the center of the flow 
cross-section are subject to higher centrifugal forces than 
those fluid particles traveling along the channel's outer 
boundary wall. As a result, the fluid in the center of the 
channel moves outwardly toward the concave boundary. As 
the process continues, the fluid particles near the boundary 
wall, move in a spanwise direction, and finally radially 
inward replacing the outwardly moving particles. These 
particles then come under the same centrifugal force and the 
process repeats itself. The resulting cyclic motion causes 


the formation of the counter rotating Taylor-Gortler 
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vortices, considered to be primarily a laminar flow 
phenomenon. 

It has been observed (Ref. 4] that the heat transfer 
rate from flow along a concave curved wall is greater than 
that for flow along a straight wall, a phenomenon attributed 
to the additional mixing provided by the secondary motion of 
the Taylor-Gortler vortices. 

There are many possible applications that could result 
from a more thorough understanding of Taylor-Gortler 
vortices and their effect on heat transfer and fluid flow 
characteristics. Two such applications could lead to 
improved heat exchanger designs and improved turbine blade 


cooling Refs. 5, 6, 7]. 


В. HISTORY 

The stability of an inviscid fluid flowing past a curved 
boundary, was first considered by Lord Rayleigh in 1916 
[Ref. 8]. By assuming that the fluid was non-viscous, he 
determined that for the motion to remain stable, its 
circulation must increase with increasing radius. С. Г. 
Taylor, in 1923, [Refs. 1, 9}, continued this study with an 
extensive analytical and experimental study of viscous 
fluids. His investigations focused on the flow between two 
cylinders, in which the inner cylinder rotated while the 
outer cylinder remained stationary. Taylor ascertained that 
such flows become unstable when the value of a dimensionless 


parameter exceeded a critical value of 41.3. 
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The parameter, known as the Taylor number is defined as: 


Та = Ве /d/Ri 


where 'd' is the width of the gap, assumed small when 
compared to 'Ri', the radius of the inner cylinder, and 'Re' 
is the Reynolds number based on the peripheral velocity of 
the inner cylinder. Taylor determined that for those cases 
in which the value of the Taylor number exceeded the 
critical value, a secondary motion developed and the Taylor 
vortices formed. Figure 2 illustrates this fluid motion. 
Instability of a similar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
first considered analytically by W. R. Dean [Ref. 10] in 
1928, for a channel formed by two concentric cylinders, 
where the radius of the inner cylinder was large in 
comparison to the spacing between the inner and outer 
cylinder walls. Dean concluded that there would be an 
initiation of flow instability, and the formation of 
Vortices, when a dimensionless parameter, the Dean number, 


exceeded a value of 36. The Dean number is defined as: 


A d 
De = Ке “БТ 


where 'd' represents the channel half-width, 'Ri' is the 


inner cylinder radius, and 'Re' is the Reynolds number based 
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on the mean velocity of the undisturbed flow and the channel 
half-width 'd'. The analytical work of Dean was later 
verified by W. H. Reid [Ref. 11], using an approximate 
Solution. 

In 1940, H. Gortler [Ref. 2] studied the stability of 
laminar boundary layer profiles on curved walls under the 
influence of small disturbances. He found that these 
disturbances were similar to the vortices studied by G. I. 
Taylor. Using approximate numerical calculations, Gortler 
concluded that the disturbances of vortices were produced 
only on the concave boundary walls and that the overall flow 
profile appeared to remain laminar in nature. These results 
were verified, with an exact solution, by G. Hammerlin, and 
reported by H. Schlichting [Ref 12 T9555. О 
Smith completed an even more extensive numerical analysis 
that further substantiated these findings [Ref. 3]. The 
results of these numerical solutions have recently been 
demonstrated with the use of hot wire anemometry, laser 
doppler systems, and flow visualizations techniques, [Refs. 
13, 14]; and in 1976, Y. Aihara [Ref. 15] conducted a non- 
linear analysis of Gortler Vortices. 

As interest began to develop concerning the effects of 
the secondary flows associated with the Taylor-Gortler 
vortices, many studies were published concerning the 


influences of these vortices on heat transfer in the laminar 
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and turbulent flow regimes. Іп 1955, Е. Ктетећ, (Кеб. 16), 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent flows. 
He concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and 
under similar turbulent conditions. 

In 1965, L. Persen [Ref. 17], considered the special 
cases of very high and very low Prandtl number fluids and 
related the increase in heat transfer from a curved wall to 
the presence of the Taylor-Gortler vortices. 

There has been only a limited amount of published 
literature dealing with the flow and heat transfer in curved 
channels of rectangular cross-section and large aspect 
ratios, where the aspect ratio is equal to the spanwise 
distance of the channel divided by the channel height. Much 
of what has been published involves the development of 
numerical approximations and solutions for heat and mass 
transfer in curved ducts of various geometries. К. Cheng 
and M. Akiyama [Ref. 18] developed a numerical solution for 
forced convection heat transfer with laminar flows in curved 
channels of rectangular cross-section, but only for small 
aspect ratios. 

In 1976, А. A. Shibani and M. N. Ozisik [Ref. 19], using 


matched asymptotic expansion techniques for a wide range of 
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Prandtl numbers, solved the heat transfer problem between 
straight parallel plates with turbulent flow, for the case 
of uniform wall temperature. 

Other experimental and analytical studies that are worth 
citing with regard to this present study are Y. Mori [Ref. 
20], who obtained results for hydrodynamically fully 
developed flows with constant wall heat flux in curved 
channels of square cross-section. W. M. Kays and E. Y. 
Leung [Ref. 21] reported solutions for turbulent flow heat 
transfer in a concentric circular tube annulus with a fully 
developed velocity profile and constant heat rate per unit 
length, for a fluid of Prandtl number 0.7. P. F. Brinich 
and R. W. Graham [Ref. 22] reported results for turbulent 
flows in a rectangular curved channel with an aspect ratio 
of 6, for the inner wall heated, the outer wall heated, and 
both walls heated. 

In 1973, at the Naval Postgraduate School, R. J. Mckee 
[Ref. 23] confirmed the presence of Taylor-Gortler vortices 
in a curved rectangular channel of aspect ratio 40. M. 
Durao [Ref. 24] modified R. J. Mckee's channel and modeled 
it as infinite parallel plates with the outer wall heated 
and inner wall adiabatic. M. Durao and J. Ballard [Refs. 
24, 25] investigated the effects of Taylor-Gortler vortices 
on the heat transfer in a straight and curved test section 


for laminar flows. В. Но11Рап, Jr. [Ref. 26] reported 
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results for laminar and transition flows. 5. Daughety 

(КеЁ., 27] reported results for turbulent flows and J. Wilson 
[Ref. 28] reported results for transition and turbulent 
flows. Each time Mckee's apparatus was used it was modified 
slightly using prior recommendations to improve on the 


accuracy of the results. 
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IL. ТЫТЕМТ ОЕ ЮУ 


The purpose of this investigation was to build an 
improved curved rectangular duct of large aspect ratio which 
could be configured with a variety of boundary conditions, 
and yield more accurate data then previously attained. By 
combining new ideas and experimental equipment with past 
experiences and recommendations of R. Mckee and others 
previously mentioned [Refs. 23, 24, 25, 26, 27, 28], the 
effect of streamwise curvature on the heat transfer rate in 
a curved rectangular channel could be studied in greater 
detail while providing a smaller uncertainty in the results. 

Some of the improvements incorporated in the new channel 
include: the ability to heat one or both walls over a wide 
range of constant heat fluxes; the ability to heat both 
walls at different constant heat fluxes; the use of more 
calibrated thermocouples which will allow greater precision 
in reading surface and air Stream temperatures, combined 
with a better arrangement which will allow the calculation 
of the local heat transfer coefficient and Nusselt number; 

a longer transition region at the outlet which will allow 
for more mixing and a more accurate outlet temperature 
measurement. The use of an outer and inner section allowed 


both wall surfaces to be specially constructed with smooth 
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transitions between wall and heater surfaces in both 
straight and curved test sections, preventing unexpected 
tripping from laminar to turbulent flow. Finally, by 
constructing the channel with materials having low thermal 
conductivities and avoiding any heat sinks along the channel 
walls, heat losses could be minimized. 

Taylor-Gortler vortices are known to enhance heat 
transfer in curved channels for transition and turbulent 
flows. This enhancement in the transfer of heat has been 
attributed to the secondary flow velocity components of the 
Taylor-Gortler vortices. These secondary components 
transport the heated fluid from the outer concave wall, 
inward toward the opposite wall of the channel, displacing 
the cooler fluid particles and causing them to move toward 
the heated concave wall. It was expected that similar 
results would be observed in this study. 

This investigation was conducted using a single channel 
with a rectangular cross-section and constant aspect ratio, 
again defining aspect ratio as the spanwise distance divided 
by the channel height. The channel incorporates a straight 
and a curved test section. Each test section consists of 
opposing walls which can be independently heated. The 
results obtained in each test section at approximately the 
same flow rates were compared in an effort to determine the 


effects of the Taylor-Gortler vortices on the heat transfer 
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rate. Also, the results of both the straight and curved 
outer sections were compared to the results of Daughety 
[Ref. 27], and Wilson [Ref. 28] to compare the new channel 
to the old one. 

The straight section results of this study were compared 
to the Dittus-Boelter correlation using hydraulic diameter 
[Ref. 29], and the results of Kays and Leung for turbulent 
and transition flow in annular passages [Ref. 21]. 

The curved section results of this study were compared 
with the results of Brinich and Graham [Ref. 22] for 
turbulent flow in a rectangular curved channel, and with 
Kreith [Ref. 16] for turbulent flows in concave and convex 


curved channels. 
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IIT. TEXPERIMENTAL: SETUP 


A. DESCRIPTION OF THE APPARATUS 

To meet the objectives of this investigation, a channel 
of rectangular cross-section and constant aspect ratio was 
constructed. The rectangular channel has four electrically 
heated wall sections, composed of a rubber heater bonded to 
a thin copper plate. All four heated wall sections are 
30.48 centimeters long and 25.4 centimeters wide, with a 
heated area of 774.2 square centimeters. The rectangular 
channel is 0.635 centimeters high and 25.4 centimeters wide 
resulting in an aspect ratio of 40. The cross-sectional 
area of the channel is 16.13 square centimeters. The wetted 
perimeter is 52.07 centimeters, and the hydraulic diameter 
is 1.239 centimeters. Figure 3 shows a cross-sectional view 
presenting both the straight and curved sections of the 
channel. 

The channel, shown in Figure 4, has a straight entrance 
region of 76.20 centimeters, which insures hydrodynamically 
fully developed flow prior to entering the straight heated 
test sections. The entrance region is followed by two 
parallel 30.48 centimeter straight heated test sections, 
which oppose each other, one in each wall. Following the 


Straight heated test sections is a 15.32 centimeter straight 
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exit region, that completes the straight test section. The 
straight test section is immediately followed by the curved 
test section which subtends a 180 degree arc, with the inner 
wall having a radius of curvature of 29.7 centimeters. The 
curved test section has a 52.63 centimeter unheated convex 
curved section prior to the curved heated sections. The two 
concentric curved heated sections, inner one convex and 
outer one concave, are 30.48 centimeters long and the outer 
heated wall subtends an arc of 58.8 degrees. The curved 
heated sections are followed by a 10.16 centimeter convex 
curved exit region which completes the curved test section. 
Finally, the channel is completed by a 92.0 centimeter 
straight exit region which allows sufficient mixing prior to 
reading the outlet temperature. The entire channel was then 
thermally insulated with a plastic foam material. 

An entrance bell constructed of Plexiglas is connected 
to the inlet of the straight entrance section of the 
channel. It was designed and manufactured according to ASME 
nozzle standards, with an elliptical curved base on a major 
axis equal to 25.4 centimeters and a minor axis of 2.54 
centimeters. A fiberglass filter was attached to the 
entrance nozzle to prevent foreign matter from entering the 
channel. An aluminum exhaust nozzle was attached to the 
exit of the channel and directed the flow from the channel 


into a two inch diameter pvc pipe. The piping contained a 
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mass flow rate measuring section which could be fitted with 
standard ASME concentric orifices. 

A total of one hundred and eighty-nine copper-constantan, 
30 gauge thermocouples are located throughout the test 
apparatus. Only ninety-nine thermocouples are read at 
anyone time, through seventy-five data channels, to record 
the desired temperatures. Each heated wall section contains 
thirty thermocouples which are read individually and then 
averaged to give the wall section surface temperature. 
Four thermocouples span the width of the channel at the 
entrance and exit to give the average entrance and exit air 
temperatures. Three groups of five thermocouples are 
connected in parallel and inserted between each layer of 
insulation at each of the heated test sections. Each group 
of five thermocouples only uses one data channel, and these 
extra thermocouples comprise the difference between the 
number of thermocouples read and the data channels 
avallable. Finally, there is one thermcouple to measure the 


air temperature at the orifice. 


B. DETAILS OF THE CONSTRUCTION 

Each of the four heated wall sections was composed of a 
copper plate laminated to a wire wound silicone rubber 
heater. The rubber heater was made by Watlow Industries. 
There are thirty copper-constantan thermocouples placed 


between the copper plate and rubber heater. Figure 5 shows 


29 


15 13 
TES кетеме 


30.48 cm 


Thermocouples 





Ш-------: 25.4 cm EU ses шј 
.635 cm ——9 


> __ 635 ст 
Outlet 
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inner plates only. 


Legend: 


Figure 5. Thermocouple Placement in Copper Plates. 
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the placement of the thermocouples on the back of the copper 
plate. The copper plate is 30.48 centimeters long, 26.72 
centimeters wide, and 0.155 centimeters thick, but only the 
middle 25.4 centimeters of copper plate width is heated and 
centered in the channel. The back of the copper plate has a 
pattern of slots milled 0.076 centimeters deep. The pattern 
of slots for the two inner plates is different from the 
pattern of slots for the two outer plates. See Figures 6, 
7, 8, and 9. The calibrated thermocouples were numbered and 
expoxied into the slots, as shown in Figures 10 and 11. The 
rubber heater has the same length and width dimensions as 
the copper plate, but only heats the middle 25.4 centimeters 
of width. This results in a heated wall area of 774.2 
square centimeters. The excess copper plate is sandwiched 
between the wall and the side spacer thereby helping to 
support the copper plate in the channel. Dow Corning 3140 
RTV was used to adhere the rubber heater to the back of the 
Copper plate, and complete the heated section. Figure 12 
and 13 represent a completed test wall. 

The channel was built on a carefully constructed wood 
frame with steel rails, that also was used to support the 
entire channel. The inner wall was composed of two 0.343 
centimeter thick and 29.85 centimeter wide pieces of 
Plexiglas laminated together, that had spaces milled out for 


the heated copper sections. The inner wall was constructed 
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Figure 6. Photograph of Milled Thermocouple Slots. 





Figure 7. Photograph of Thermocouple Placement Holes. 
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Figure 8. Photograph of Both Inner Plates Showing 
Thermocouple Slot Pattern. 





Figure 9. Photograph of Both Outer Plates Showing 
Thermocouple Slot Pattern. 
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Figüre T2: 


Photograph of Completed Heated Test Section. 
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by heating and forming the Plexiglas around the wood frame. 
The Plexiglas wall was secured to the frame by screws tapped 
into the steel rails. The heated copper sections were then 
fitted into the machined spaces and two 0.635 centimeter 
thick and 2.223 centimeter wide Plexiglas spacers, that 
Served as the sides of the channel and support for the 
copper plates, were attached by screws to the inner wall and 
Steel rails. The two spacers were 25.4 centimeters apart 
and formed the width of the channel. Additionally, five 
flat head screws spanning the width of the channel at the 
joints where the copper and Plexiglas met were used to 
secure the copper plates and maintain an even joint. The 
Seams were then worked to a smooth transition by applying 
"Bondo" to the joints and sanding. See Figures 14, 15, 16, 
and 17. 

The outer wall was constructed in a similar manner, but 
using 0.318 centimeter thick and 29.85 centimeter wide 
pieces of Lexan.  Lexan was used on the outer wall because 
of its flexibility and ease of bending without heat, which 
would be required when attaching the outer wall to the inner 
wall and frame. Also, the outer wall was formed in an outer 
wood frame so that the transition joints of the outer wall 
could be worked smooth like the inner wall. See Figures 18 
and 19. Both walls were constructed so that no seams or 


butt joints existed on the inside of the channel, except 
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Figure 14. Photograph of Channel Frame. 





Figure 15. Photograph of Bottom View of Inner Straight 
Test Section. 
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Figure 16. Photograph of Inner Straight Wall at 
Copper-Plexiglas Joint. 





Figure 17. Photograph of Inner Curved Wall at 
Copper-Plexiglas Joint. 
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Photograph of Outer Frame and Lexan Wall 


Showing Straight Test Section. 


Figure le. 
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Photograph of Outer Curved Test Section. 


Figure 19. 
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where the heated copper section joined the Plexiglas and 
Lexan walls, and all eight of those joints had a smooth 
transition. 

‘Once both walls were completed, the outer wall and 
support frame was positioned onto the inner wall and main 
frame. The outer wall was aligned onto the inner wall and 
the outer wall support frame removed. The outer wall was 
then bolted to the inner wall and frame with steel bars 
running the length of the channel used as a spacer to insure 
even compression. Once this was completed the inlet and 
outlet nozzles were attached to the channel. 

The sides of the channel and all connections between the 
channel, the pvc piping, and the orifice were sealed with 
Macklanburg-Duncan Silicone Rubber Sealant Caulk to ensure 
there was no leakage of air into the channel or the piping 
which could effect the recorded temperatures or mass flow 
rate. The entire channel was thermally insulated with three 
layers of 1/2 inch Armstrong Armaflex 22 Sheet Insulation, 

a flexible foamed plastic material. By positioning 
thermocouples between these insulation layers, the heat lost 
through the insulation could be computed at each test 
section. The heat loss through insulation surrounding the 
rest of the channel due to the heated air flowing through 
the channel was a small fraction of the total power supplied 
to the heater. The insulation was held in place by the use 


of silver duct tape and velcro strips attached to the frame. 
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C INSTRUMENTATION 

The data acquisition system used for this experiment was 
a Hewlett Packard 3054A Automatic Data Acquisition/Control 
System consisting of a 3456A Digital Voltmeter, a 3497A Data 
Acquisition/Control Unit, and four input panels capable of 
reading any mix of eight voltages or copper-constantan 
thermocouples. Also used in conjunction with this data 
acquisition system, were a Hewlett Packard 9826 Computer and 
2671A Printer. The program used to reduce the data is 
listed in Appendix A. 

The copper-constantan thermocouples were calibrated 
against a platinum resistance thermometer using a ROSEMONT 
Commutating Bridge model 920A, and a ROSEMONT Constant 
Temperature Bath model 913A. All thermocouples were 
calibrated between 20 and 85 degrees Celsius, and a second 
order polynomial was curve fitted to correct each 
thermocouple. 

The mass flow rate in the channel was measured with a 
thin plate orifice measuring device which was placed in the 
2 inch pvc piping between the channel and the compressor. 

In this study an orifice with a diameter of 2.7318 

centimeters was used to measure hydraulic Reynold's numbers 
above 5000. The orifice device was installed in accordance 
with ASME Power Test Code [Ref. 30], with 3/8 inch pressure 


taps placed on either side of the orifice at 1 апа 1/2 
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diameters, and one thermocouple inserted in the two inch pvc 
piping down stream of the orifice to record the temperature 
of the air flowing through the orifice. Pressure measure- 
ments were taken with Meriam vertical manometers, one whose 
fluid was water with a 0 to 60 inch range, another using 
water with а 0 to 30 inch range, and a third whose fluid 

was mercury, calibrated to read inches of water with a range 
of 0 to 415 inches. Atmospheric pressure was measured with 
a Princo precision barometer. The pressure readings were 
used to calculate the mass flow rate of the working fluid as 
outlined in [Refs. 30, 31]. 

Power was supplied to the heated test section by using 
two Lambda Model LK345A FM 60 volt power supplies in series. 
The power supplies are rated at 6.0 amps. maximum at 40 
degrees Celsius. Since the electrical resistance of the 
rubber heater is not constant, but varies slowly with 
temperature, a precision resistor was calibrated and 
connected in series with the heater to allow the calculation 
of the instantaneous power being supplied. Both the outer 
and inner heated test sections have their own circuitry 
and precision resistor. The value of the precision 
resistor for the outer test sections is 2.0173 ohms, and 
the value for the inner test sections is 2.0084 ohms. See 
Figure 20 for diagram of circuit used to measure 


instantaneous power. 
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Figure 20. Power Measurement ЕТЕ 
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The working fluid was air at room temperature, which was 
drawn through the channel by an electrically driven Spencer 
Turbo Compressor, rated at 30 horsepower at 3500 rpm, and 
550 cubic feet per minute at 70 degrees fahrenheit and one 
atmosphere. Photographs of the channel and associated test 


equipment are shown in Figures 21, 22, 23, and 24. 
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Photograph of Data Acquisition System. 
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Photograph of Flow Measuring Apparatus. 


Figure 22. 
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Photograph of Spencer Turbo Compressor. 
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Photograph of Test Channel. 


Figure 24. 
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ІУ. EXPERIMENTAL PROCEDURES 


The experimental procedures followed were the same for 
the straight and curved test sections. The procedure for 
Starting up the experimental apparatus is as follows: 

(1) turn on the computer and data acquisition system 
(2) close globe valve at channel exit 

(3) close valve to manometer reading Pl 

(4) set compressor blast gate to 1/8 open 

(5) start turbo compressor 

(6) open intake valve on compressor 1/4 Ко 1 turn 
(7) open globe valve at channel 


(8) open manometer valve slowly to ensure Pl does not 
exceed manometer scale 


(9) turn on and adjust power supply to the heater 


(10) use intake valve on compressor to control mass flow 
rate through the channel 


It was necessary to determine the approximate voltage 
setting required to bring the heated wall to approximately 
50 degrees Celsius. This heated wall temperature would 
ensure a large temperature difference between the inlet and 
outlet temperatures, and the heated and unheated wall 
temperatures, while not damaging the apparatus. The larger 
the temperature difference the smaller the uncertainty in 


the results. The required power setting was found by 


51 


monitoring the heated plate temperature locally from the 
data acquisition system while the voltage to the heater was 
adjusted. 

To compute the instantaneous power supplied to the 


heater the following relationship was used: 


Уы Урв 


Q = 
P RoR 


The precision resistor voltage and the heater voltage were 
read by the data acquisition system and the resistance of 
the precision resistor was a known constant. А diagram of 
the circuit used was shown in Figure 16. 

The mass flow rate was monitored by the two manometers 
at the orifice. For each run, the atmospheric pressure, the 
pressure difference across the orifice, and the pressure 
upstream of the orifice were entered manually into the 
computer for the mass flow rate calculation. All other data 
was acquired by the data acquisition system directly from 
the experimental apparatus. 

Once the apparatus was started properly, sufficient time 
was required to reach steady state prior to gathering data. 
Preliminary runs were performed to determine the time 
required for the test apparatus to reach steady state. For 
this experiment at least two hours were allowed for a cold 


test apparatus to reach steady state, and at least one hour 
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if the apparatus was already running. After waiting the 
required time, data was taken at ten minute intervals to 
ensure that a steady state condition had been achieved. The 
criteria for steady state was based on three variables: 

(1) the difference between inlet and outlet temperature 

(2) the heated test section wall temperature 

(3) the unaccounted heat loss in the system 
When the first two of these quantities varied by less than 
two percent over a ten minute interval and the third was 
less than ten percent of the input power, it was considered 
that a steady state condition had been achieved. 

Experiments were conducted at specific mass flow rates 
of air corresponding to Reynolds numbers between 7000 and 
22000. The Reynolds number was based on the hydraulic 
diameter of the channel. For each mass flow rate, tempera- 
ture and voltage data were taken automatically by the data 
acquisition system and immediately reduced and printed. The 
variables that were measured and used in this investigation 
were: 

(1) the temperature of the air entering the channel Өв.) 
В temperature of the air leaving the channel (Tout? 


(3) the temperature of the heated wall in the test 


section (156 ог Te 


(4) the temperature of the unheated wall in the test 


Ест оп (Т ог а) 


(5) the temperature between each layer of insulation at 


the test Section Chee) 
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(6) 
(7) 
(8) 
(9) 


(10) 


(11) 


the 
the 
the 
the 


the 


temperature of the air at the orifice ee) 
voltage across the precision resistor (Урд) 
voltage across the rubber heater (Уң) 
pressure upstream of the orifice (P1) 


difference in the upstream and downstream 


pressure across the orifice (AP) 


the atmospheric pressure (Р, тм) 


The results obtained from the straight test sections 


Served aS a baseline for the comparison of the curved 


section results. An effort was made to match the mass flow 


rates for each test section as closely as possible. 
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V. PRESENTATION OF DATA 


A. DATA REDUCTION 

The constant heat flux surface for the heated wall of 
each test section was provided by the uniform electrical 
resistivity of the rubber heater. The insulated unheated 
wall was considered adiabatic since the heat losses through 
that wall were negligible. The straight portion of the 
channel upstream of the straight heated test section was of 
sufficient length to ensure that the air flow was hydro- 
dynamically fully developed for the flow velocities of this 
experiment. The straight portion of the channel downstream 
of the curved section proved to be of sufficient length to 
ensure that the air flow exiting the channel was thoroughly 


) 


mixed and that the measured outlet air temperature (Tawe 
was an average bulk temperature, as long as the hydraulic 
Reynolds number of the channel was less than 17000. 

The analysis of the data required several quantities to 


be defined. The heat convected to the air was calculated 


using the equation: 


ЕЕ mc 


air pair (T - Tin? 


out 1n 


where Cpair is the specific heat of the air at constant 


pressure and 'm' is the mass flow rate of the air. 
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The average heat transfer coefficient between the heated 
wall and the flow of air in the channel was defined by the 


equation: 


F =. Qair 
АРІ, АТ 
і ! а 2 ' і > 
where D 15 defined above, АРІ, 1s the area of the 


heated copper plate in the test section, and '^AT' is the 
difference between the arithmetic average heated wall 
temperature CEF) and the average bulk air temperature 
(Tpik)* The average bulk air temperature is defined as the 
arithmetic mean of the entrance and exit air temperatures 


CT T ды 


in^ gut 


The average Nusselt number was then calculated using: 


air 


In this equation "Dig. is the hydraulic diameter and Ки 
is the thermal conductivity of the air. 
The hydraulic Reynolds number was calculated for each 
test run as follows: 
m D 
Re = He 


hd А р. 
C alr 


where again 'm' and "Бра! are the mass flow rate and 


hydraulic diameter of the channel, Ac! is the cross- 
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sectional area of the channel, and E is the dynamic 


viscosity of the air. 


For the curved sections the Dean number was defined as: 





цо”, С 
e A u. Ri 2 Е1 
c air 
where "Кез! is the Reynolds number based on the channel 


height, p is the channel height, and 'Ri' is the radius 
of curvature of the inner convex wall. 

A Sketch of the control volume and a set of sample 
calculations for one test run of the curved section are 
given in Appendix B. А sample calculation for the 


uncertainty analysis is given in Appendix C. 


B.  RESUCTS 

The data obtained from each experimental run was reduced 
by the computer program listed in Appendix A, utilizing the 
expressions described in Part A above. The major parameters 
resulting from this process are shown in Tables I through 
IV. Tables I and II contain the straight sections results 
and Tables III and IV contain the results of the curved 
sections. A plot of the average Nusselt number versus 
hydraulic Reynolds number is given in Figures 25, 26, 27, 
28, and 29 for the comparison of the results. "Error" bands 


have been indicated as a result of the uncertainty analysis. 


БИ 


Ке а 
7709 

11953 

15646 


21455 


и. 
(Я) 


198-50 
135.26 
157592 


дані З 


TABLE I 


(W/mC ) 


41.99 
ои 57 
74.43 


87.42 


АТ 
(С) 
285.2 
25.49 
24.45 


27.36 


SUMMARY OF STRAIGHT OUTER WALL TEST RESULTS 


hd 


20.07 


ле >т! 


32465 


41.95 


Reynolds number and Nusselt number are based on 


hydraulic diameter. 


Results for Reynolds number of 


21455 based on orifice air temperature instead of 
outlet air temperature. 
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7628 


12169 
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ШЕВ 256 
159.84 


204.42 


TABBE 11 


h 


(W/mC) 


41.52 


59. 50 


70.24 


rece 


AT 
(C) 
28239 
27.42 
26.37 


27.84 


SUMMARY OF STRAIGHT INNER WALL TEST RESULTS 


Мара 
19.85 
28.50 
ЗЭ?Л 


393 1 


Reynolds number and Nusselt number are based on 


hydraulic diameter. 


Results for Reynolds number of 


21312 based on orifice air temperature instead of 


outlet air temperature. 
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Вена 

7782 
11726 
15813 


21246 


De 


206 
ЗОРЧ 
419 


563 


8}: 


115.88 
157.45 
158. 17 


231.39 


TABLE III 


h 


(W/mC) 


52.43 
68.16 
Ва 


93.02 


(С) 


25:51 
PITT 
22.34 


27.20 


SUMMARY OF CURVED CONCAVE WALL TEST RESULTS 


Мана 
25.00 
32753 
39.78 


44.52 


Reynolds number and Nusselt number are based on 


hydraulic diameter. 


Results for Reynolds number of 


21246 based on orifice air temperature instead of 


outlet air temperature. 
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TABLE IV 


SUMMARY OF CURVED CONVEX WALL TEST RESULTS 


Re De @ һ АТ Nu 
hd h 
CM (W/mC) (С) З 
7757 206 82.37 38: 205207 18.37 
12021 318 ДЕ 50 57.62 23.92 28 . 065 
15908 421 122524 68.57 240. 34 32.94 
21603 5 72 164.08 79.28 23243 33.11 


Reynolds number and Nusselt number are based on 
hydraulic diameter. Results for Reynolds number of 
21603 based on orifice air temperature instead of 
outlet air temperature. 
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As would be expected the results indicate an increase in 
the rate of heat transfer with increasing Reynolds number 
for both the straight and curved test sections. In 
addition, the heat transfer rate was highest in the concave 
curved section for each Reynolds number investigated, and 
lowest in the convex curved section. 

In the low end of the turbulent regime, at a hydraulic 
Reynolds number of 7780, the rate of heat transfer was about 
twenty-two percent higher in the concave curved section than 
the outer straight section. As the hydraulic Reynolds 
number continued to increase, between 11700 and 15800, the 
rate of heat transfer was approximately eleven percent 
higher between the concave curved section and the outer 
Straight section. At a hydraulic Reynolds number of 21000 
the difference in heat transfer rate decreased to about 
семе" ретсепь, 

The results obtained at a hydraulic Reynolds number of 
21000 have a much larger uncertainty, because the orifice 
air temperature was higher than the average outlet air 
temperature. This signifies that the air was not fully 
mixed at the outlet, because no heat is added to the air 
after the outlet air temperature is read, yet the air 
temperature increased. For hydraulic Reynolds numbers 
greater than about 17000 the channel outlet section is not 


long enough to allow the air to mix properly. 
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The values obtained for the hydraulic Reynolds number of 
approximately 21000 are based on the orifice air temperature 
vice the outlet air temperature. This resulted in a Nusselt 
number higher than obtained by using the the outlet air 
temperature, but still lower than the actual Nusselt number, 
because of the small heat loss between the outlet nozzle and 
the thermocouple at the orifice. This heat loss cannot be 
accounted for and the temperature used as an outlet air 
temperature was still lower than the actual outlet air 
temperature. 

In comparing the concave curved wall to the convex 
curved wall the rate of heat transfer was about thirty-six 
percent higher in the concave section for a hydraulic 
Reynolds number of 7780. For hydraulic Reynolds numbers 
between 11700 and 15800 the rate of heat transfer was twenty 
percent higher in the concave section, and at a hydraulic 
Reynolds number of 21200 the rate of heat transfer was 


Seventeen percent higher. 
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VI. DISCUSSION AND CONCLUSIONS 


J. Wilson [Ref. 28], in his study, reported that the 
actual presence of Taylor-Gortler vortices were verified by 
the observation of the liquid crystal isotherm distribution 
in the curved test section. Не further reported, that at 
low values of hydraulic Reynolds number and through most of 
the transition regime, liquid crystal bands were oriented in 
the streamwise direction. As the Reynolds number increased 
turbulent forces began to dominate in the curved section and 
less evidence could be detected for the existence of the 
vortices. 

The high aspect ratio of the channel provided the basis 
for the assumption that the experimental apparatus should 
behave like infinite parallel plates. Holihan's data [Ref. 
26] tended to verify this assumption. His experimental data 
in the laminar flow region, approached the theoretical limit 
for average Nusselt number for parallel plates with one wall 
heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 32]. Based on the assumptions mentioned 
above, comparisons were made with other analytical and 
experimental results that were of the same or similar 
problems. 

A comparison of the experimental results of this study 


and the experimental results of Daughety [Ref. 27] and 
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Wilson [Ref. 28] are shown in Figure 30 for the curved 
section, and Figure 31 for the straight section. The 
channel aspect ratio was 40 for each of these studies and 
the experimental procedures were similar, with the exception 
of a slight modification in Wilson's procedures. Wilson 
[Ref. 28] made a correction to determine the temperature 
drop across the liquid crystal sheet covering the heated 
wall. By applying this temperature drop to the measured 
temperature, the true temperature of the heated wall would 
be known. After reviewing Wilson's data and experimental 
procedure it was detected that Wilson was over estimating 
the temperature drop across the liquid crystal by 
approximately five to seven percent, and this lead to an 
over estimation of T of ten to twelve percent. This 
resulted in an over estimation of the average heat transfer 
coefficient and the average Nusselt number of between ten 
and twelve percent. When Wilson's data is lowered by ten 
percent, both his data and Daughety's data fall within the 
uncertainty bands of the present study as shown in Figures 
32 and 33. Differences were attributed to the inability to 
acccount for all of the heat transfer processes and the 
limitations inherent in any experimental work. 

Ballard [Ref. 25] reported an increase in the rate of 
heat transfer between a concave wall and a straight outer 
wall of eleven percent for laminar flows.  Holihan [Ref. 26] 


reported an increase of fifteen percent for laminar flows 
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and thirty percent for transition flows.  Daughety [Ref. 
27], reported an increase of twenty percent in the heat 
transfer rate for turbulent flows. Wilson [Ref. 28], 
reported an increase of forty-five percent for laminar 
flows, twenty percent for transition flows, thirteen percent 
for hydraulic Reynolds numbers between 5400 and 12000, and 
about twenty-six percent at a hydraulic Reynolds number of 
24000. This compares very favorably with the present 
results, except for the data at a hydraulic Reynolds number 
of 21000. Мо conclusion can be made about the present 
results at a hydraulic Reynolds number of 21000 because of 
the high uncertainty of the data as previously mentioned. 
This area needs further experimentation before accurate 
conclusions can be made about the heat transfer rate at high 
hydraulic Reynolds numbers. 

Kreith [Ref. 16] reported an increase in the rate of 
heat transfer along a concave wall as compared to a convex 
wall of twenty-five to sixty percent for hydraulic Reynolds 
numbers between oS to іне This also compares favorably 
with the present results which indicate an increase of 
seventeen to thirty-eight percent when comparing concave 
wall to convex wall. 

A comparison of analytical and experimental results with 
those of the present study for the straight test section is 


shown graphically in Figure 34. The Dittus-Boelter equation 
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[Ref. 29] for heat transfer in a straight tube for constant 


wall temperature is plotted from: 


Ми = 0.02Rne9- 8 


For a Prandtl number equal to 0.71. Additionally, the 
results of Kays and Leung [Ref. 21] for heat transfer in 
annular passages are shown for r* - 1.0, which equates to 
parallel plates. When examined these results compare 
favorably with the present data. 

The curved section data is plotted in Figure 35, and is 
compared to the experimental work of Brinich and Graham 
[Ref. 22]. The accuracy of the data points used from 
Brinich and Graham are subject to some error, in that the 
actual values are not given in their study and had to be 
taken from a plot of Stanton number versus Reynolds number. 

After carefully examining the past and present data the 
trends and results are consistent. Heat transfer rates 
increase with increasing Reynolds numbers. Also, the 
сопсауе curved section enhances the rate of heat transfer 
over the convex curved section and the straight test 
sections for laminar, transition, and turbulent flows. Тһе 
presence of Taylor-Gortler vortices in the laminar and 
transition regime seems to be the major contributing factor 
to this enhancement. Whether the Taylor-Gortler vortices 
continue to exist in the turbulent flow regime remains to be 


Seen. 
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VII. RECOMMENDATIONS 


The amount of experimental work remaining in the study 
of heat transfer in curved channels is still significant. 
The present experimental apparatus has been tested in the 
turbulent regime, and the test data proves that accurate 
results can be obtained for hydraulic Reynolds numbers below 
17000, slightly higher when heating the straight test 
section. Additional experiments into the laminar, 
transition, and turbulent region should be conducted to 
better understand the heat transfer processes occurring in 
the straight and curved sections. By remaining below a 
hydraulic Reynolds number of 17000 accurate data can be 
obtained and significant conclusions can be drawn in the 
comparison of the curved sections to the straight section, 
and the concave curved wall to the convex curved wall. 

To improve the accuracy of the data and results some 
Slight modifications to the present channel should be done. 
First, the straight outlet section of the channel should be 
lengthened to ensure proper mixing of the heated outlet air. 
Second, some sort of honey combing, screening, or other type 
of device should be placed in the outlet section to enhance 
the mixing of the heated outlet air. Finally, extra 


thermocouples should be placed at the outlet and read more 
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then once during a data run to decrease the uncertainty of 
the outlet air temperature E or! and the uncertainty of the 
difference between the inlet and outlet air temperature. By 
applying these recommendations the uncertainty of the 
Nusselt number can be reduced to four percent or less. 

Further experiments at hydraulic Reynolds numbers 
greater than 17000 also need to be examined. This can be 
accomplished if the above recommendations are taken along 
with one other. The channel walls need outer support to 
prevent the center of the channel from collapsing at the 
high hydraulic Reynolds numbers. Slight collapsing was 
noticed at hydraulic Reynolds numbers greater than 20000. 
This collapsing was caused by the laminating of two 0.3175 
centimeter thick walls which were not as rigid as the single 
0.635 centimeter thick walls used previously. The actual 
hydraulic Reynolds numbers that can be attained and still 
yield accurate results will be limited by the power of the 
compressor and the rigidity of the channel's walls. 

Finally, once a strong data base can be obtained for 
each test section further experimentation can be conducted 
with different boundary conditions. For example, both 
curved walls heated at the same constant heat flux. Once 
this can be done a true comparison of the heat transfer rate 
of the total curved test section can be compared to the 
total straight test section, i.e., both curved walls versus 


both straight walls. 
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APPENDIX A: COMPUTER PROGRAM 


! FILE NAME: DRP2A 
LAST REVISION DATE: 21 OCTOBER 1985 
! AUTHOR: LT. GEORGE С. GALYO 


PROGRAM FOR GATHERING AND REDUCING DATA WHEN HEATING ONE COPPER PLATE, 
WITH 10 MINUTES BETWEEN RUNS. PROGRAM IS MODIFIED TO ACCOUNT FOR 
BROKEN THERMOCOUPLE CHANNELS 5, 16, 21, AND 26 WHICH EFFECT OUTER 
еер HEATER, AND THERMOCOUPLE CHANNEL 47 WHICH EFFECTS INNER CURVED 
HEATER. 


COM /Co/ D(7),Aa(76,2) 
DIM Emf(79),T(79) 


! CORRELATION FACTORS FOR EMF TO DEGREES CELSIUS FOR CU-CO T/C. 


DATA 0.10086091,25727.94369,-767345.8295.78025595.81 

DATA -9247486589,6.97688E+11,-2.66192E+13,3.94078E+14 

READ р(+) 

PRINTER IS 701 

BEEP 

INPUT "ENTER MONTH, DAY, AND TIME (MM:DD:HH:MM:SS) *,TimeS 

OUTPUT ?709;"TD";TimeS$ 

BEEP 

INPUT “ENTER INPUT MODE (0=3054A,1=FILE)*,Im 

IF Im=1 THEN 

BEEP 

INPUT "ENTER NAME OF EXISTING DATA FILE" ,Oldfiles 

PRINT USING *20X,""These Results Are From Data File: "".108":01Д4611е5 
PRINT 

PRINT 

ELSE 

BEEP 

INPUT “ENTER NAME ОР МЕН DATA FILE?*,Neuf1le$ 

PRINT USING "20X,""Data For This Run Is Stored In Data File: °*,10A°:Newf 


PRINT 


PRINT 
END IF 


! ENTER THE CHANNEL CONFIGURATION. 


PSI I см NEP 


RINT USING *9X,""THE FOLLOHING CHANNEL CONFIGURATION HAS SELECTED:"*"* 
I 

BEEP 

INPUT “SELECT HEATER TYPE (O=STRAIGHT, 1=CURVED)".Itype 

IF Itype-0 THEN PRINT USING "12X,""Heatinag STRAIGHT Test Section.""" 
ae THEN PRINT USING "12X,""Heatina CURVED Test Section.""" 


INPUT "SELECT HEATER POSITION (O=OUTER, 1=INNER)".Ipos 

IF Ipos=0 THEN PRINT USING *12X,""Heating OUTER Plate."^" 
is Ipos-1 THEN PRINT USING "12X,""Heating INNER Plate.""" 
RINT 


! ENTERS FILES WITH T/C CALIBRATION COEFFICIENTS BASED ON HEATER TYPE. 


IF Itype-0 THEN ASSIGN @Filec TO “STRA_COEF" 
IF Itype=1 THEN ASSIGN (Filec TO "CURV, COEF" 
FOR I=0 TO 76 

ENTER @Filec:Aa(I.0),Aa(I.1),Aa(I,2) 

NEXT I 
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1590 
1600! 
1610 
1620! 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740! 
1750 
1760! 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890! 
1900 
1910! 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
99) 

2000 
dd) 

2010 
еш” ы 
2020 
cj 
2030 
2040 
2050 
2060! 
2070 
2080! 
2090 
2100 
2110 
2120 
2130 
2140 


ASSIGN @Filec TO * 
| INSERT DATA STORAGE DISK. 


BEEP 

INPUT "CHANGE DISK AND HIT ENTER",Ok 
IF Im=1 THEN ASSIGN @File TO Oldfile$ 
IF Im=0 THEN 

CREATE BDAT Newfile$,40 

ASSIGN @File TO Newfile$s 

END IF 

ENTER 709; TimeS 

CLEAR 709 

PRINT USING *20X,""Month, Day, and Time: "".14A";TimeS$ 
PRINT 


! FOR DATA RUN READ AND STORE ALL RAW EMF VALUES. 


IF Im=0 THEN 

OUTPUT 709;"AR AFOO AL79" 
OUTPUT 722:" FT RI T1 Z1 FL1" 
FOR I=0 TO 79 

OUTPUT 709;"AS" 

ENTER 722;Emf CID 

NEXT I 

OUTPUT @FilesEmf(+*) 

ELSE 

ENTER @FilesEmf(*) 

END IF 

OUTPUT 709;"TD" 


! CONVERT AND CALIBRATE ALL T/C READINGS TO DEGREES CELSIUS. 


hl USING "9X,""THE FOLLOWING DATA WERE RECORDED: """ 
FOR I=0 TO 79 

IF I-38 OR I-39 OR I-77 OR I-78 THEN 

IF (Ipos=0 AND I=38) OR (Ipos=1 AND I=77) THEN 


ladd=0 
IF Ipos*1 THEN Iadd-39 
PRINT USING "12X,""Heater Voltage (Vh) = "" 20.30,"" (У)"““гЕтг (ам а 


PRINT USING "*12X,""Resistor Voltage (Vr) м-”9,20:30,"" (9УУ7"”"” ЕпВМГа 
IF Ipos=0 THEN PRINT USING "12X,""Precision Resistor (Rpr) = 2.01729 (Ohm 
IF Ipos=1 THEN PRINT USING "12X,""Precision Resistor (Rpr) = 2.00839 (Ohm 


PRINT 
END IF 
ELSE 


! IF T/C OUT OF RANGE PRINT A QUE. 


Aemf=Emf(I) 

IF Aemf<3.E-4 OR Aemf>3.E-3 THEN 
1(1)=999.99 

ELSE 

CALL TvsvCEm£ CI), T£) 

IF I=79 THEN 
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T(I)=Tt 

GOTO 2210 

END IF 
Т(І)«РМТсогг (ТІ, 1) 
END IF 

END IF 

NEXT I 


2230 ! PRINT ALL T/C READINGS. 


2150 
2160 
2170 
2180 


2310 


РАТМТ 
PRINT 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


15). 1(17) 
2320 PRINT 
DIM" 


2930 РКІМТ 
2), Т (23) 
2340 PRINT 
EDS" 

2350 PRINT 
28) .T1(29) 


USING 
USING 


USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


ТЕХ UUTERSISEBDES- ^ 


*14X, ""T/C Number: 
"14X, ""Temp (C): 
“14X,""T/C Number: 
"14X, ""Temp (С): 
*14X, "*T/C Number: 
"Q1l4X,""Temp (C): 
"14X, ""T/C Number: 
"14Х,""Тетр (С): 
"14Х,""Т/С Мипбег: 
"14Х,""Тепр (С): 


1 2 3 4 5 
С СЗО ООА ТСО, ТСТ СЗУ ИИ, т 
7 8 9 10 11 
Vol SD e0D yan). 106), 107), FCS), 1. soem), 
12 14 |25; 16 17 
99,6 ВОНЕХ ^; ТО, ГОО, ТО. Т( 
19 20 21 22 23 
шаг с(2020002?285)2:1(618),1019822:21620),120,11 
25 25 27 28 29 
АЕ БІЗ С ТИС БУТСО СИЕТ. 


2360 
2370! 
2380 
2390! 
2400 
2410 
2420 
2430 
2440 
2450 
2460! 
2470 
АТЕК. 
2480! 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2550 
2600 
2610 
2620 


PRINT 


! CALCULATE AVERAGE TEMPERATURE OF OUTER HEATED WALL. 


Sum=0. 

IF Itype-0 THEN 
РОК 1-0 TO 29 
Sum*Sum*T(CI) 
NEXT I 
Two=Sum/30 


! CORRECTION FOR BROKEN T/C CHANNELS 5, 16, 21, AND 26 IN OUTER CURVED HE 


ELSE 

РОК 1-0 Т0 29 

[Е 1-5 ОКВ 1-16 ОК 1-21 ОВ I-26 THEN 2530 
Sum*Sum*T(CI) 

NEXT I 

Tuo*Sum/26 

END IF 

PRINT USING "12X,""OUTER INSULATION: "** 

PRINT USING "14X,"*T/C Number: 31 (Tinsol) 


PRINT USING "14X,""Temp (C): 
Tinso1*T(30) 

Tinso2=1(31) 
Тїп503:1(32) 
PRINT 


32 (Tinso2) 33 (110503)“ 
ЕО ВВ. Ху (ЗОВ). TC 328 
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2630 
2640 
2650 
2660 
2670 
2680 
2690! 
2700 
2710! 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
4Б""" 
2790 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


USING 


USING 
USING 
USING 


“12X. "ORIFICE TIENE Тот?) ме 27,30.00,77 


(C)" 242093) 


"12X, " "INEETSTEMPERATURE : * * * 


“14X,°"T/C Number: 
"14Х,""Тепр (С): 


35 36 37 28777 
=" 4(3D.DD,2X)"; T0347, 1035), 7 Cam, 11379 


PRINT 
! CALCULATE AVERAGE INLET TEMPERATURE. 


Sum=0. 
FOR I*34 TO 3? 
Sum*Sum*T(I) 
NEXT I 
Tin=Sum/4 
PRINT USING 
PRINT USING 


PRINT USING 


*12X, " INNERSBEBRTE:"*" 
"14Х,""Т/С Мипрег: 


"14Х, " "Тетр (С): 


41 42 43 44 45 


"5,БС30:00.2Х)":1640) 16414 (042) , Tey), T 


44),T(C45) 


2810 


PRINT 
PRINT 


USING 
USING 


“14Х,““Т/С МитБег: 47 48 49 50 9i 


"14Х,""Тетр (С): "7,5(30.00,2Х)7;Т1(46),1(47),Т1(48),Т(49),Т( 


50) , ШЕШ) 


2820 


2830 


PRINT 
PRINT 


USING 
USING 


*14X, "*T/C Number: 
"TAX." "Tempus (0): 


53 54 55 56 БА 
“*,6(3D.DD,2X) 01 (52). Т0580:, ТОБА), ТЕТ 


S6) , TOES 


2840 


2850 


PRINT 
PRINT 


USING 
USING 


"14Х,""Т/С Мипрег: 
"14X, ""Temp (C): 


59 60 61 62 63 
**,O6CD.DD,2X)* 2 T0580), TC590, ТОБ, ТОЕТ 


62),1(63) 


2860 
ZU" 
28/0 


PRINT USING 


USING 


“14Х,““Т/С Number: 65 66 67 68 баце 


PRINT "”,6(30-00,2Х)5:764),1(65),Т665),ТСе С 


"14Х,""Тетр (С): 


68), Т(69) 


2880 
2830! 
2900 
2810! 
2320 
2930 
2940 
2950 
2960 
2970 
2980! 
2990 
3000! 
3010 
3020 
3030 
3040 


PRINT 
! CALCULATE AVERAGE TEMPERATURE OF INNER HEATED WALL. 


Sum=0. 

IF Itype=0 THEN 
FOR I-40 TO 69 
Sum=Sum+T (1) 
NEXT I 
Tw1=Sum/30 


! CORRECTION FOR BROKEN T/C CHANNEL 47 IN INNER CURVED HEATER. 


ELSE 

FOR I-40 TO 69 

IF I*47 THEN 3050 
м 


3050 МЕХТ 


3060 
3070 
3080 
3090 
3100 
eus 


Twi=Sum/29 
END IF 


PRINT USING "12X,""INNER INSULATION: *** 
PRINT USING *14X,*"T/C Number: (13: Pins TT) 


PRINT USING "14X.""Temp (C): 
Тіпсі1-Т(70) 


А 23 (Тата 3)" 
""ЗОЗу. Ир, В ИО ПА Се 
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Tinsi2-T(71) 


Tinsi3=T(72) 

PRINT 

PRINT USING *12X,"**OUTLET TEMPERATURE: "*"* 

PRINT USING *14X,"*T/C Number: 74 75 76 пио“ 

PRINT USING "14Х,""Тепр (С): 00:2)" TIS) ТОР 1025). ТЕЛЬ 


PRINT USING "&,»" 

! CALCULATE AVERAGE OUTLET TEMPERATURE. 
Sum=0. 

FOR I=73 TO 76 

Sum=Sum+T (I) 


NEXT I 
Tout*Sum/4 


! CALCULATE DIFFERENCE BETWEEN CHANNEL INLET AND OUTLET TEMPERATURES. 
Tdiff\*Tout-Tin 

! CALCULATE CHANNEL BULK TEMPERATURE. 

Tblk=(TintTout) /2 

! PRINT CALCULATED TEMPERATURES. 


Por USING "SX,**THE FOLLOWING TEMPERATURES МЕКЕ САСО АТЕр: ** * 
Р 


LR ЕК по “12Х, ""Ауегаае Оифег Hall Temperature (Tuo) = =“ ар. 

‚1шО 

к ы “12X,""Average Inner Wall Temperature (Twi) ео). 

улшаг wi { | 

PRINT USING "12X,"*Average Outlet Temperature (Tout) « SED. 

(С) ЕШ ОНЕ 

RU Be ME “12X,""Average Inlet Temperature (Tin) = 7%). 
Тіп 

о Channel Inlet and Outlet Temp Diff (Tdiff) = "SED. 
> Idi 

PRINT USING *12X,""Average Bulk Temperature (Tblk) = “Жр. 


(Cy 27>, [blk 


! CALCULATE DELTA T (Tdel) FOR USE IN CALCULATING H LATER. 


Tdelo=Two-Tblk 
Tdeli=Twi-Tblk 
IF Ipos=0 THEN Tdel=Tdelo 
IF Ipos-1 THEN Tdel-Tdeli 


PRINT USING "12X,""Mean Temperature Difference (Tdel) = SED. 


(Clean Тае] 
! CALCULATE DIFFERENCE BETWEEN OUTER AND INNER WALL. 
Tuwdiff=ABS (Two-Twi ) 


PRINT USING "12X,""*Outer and Inner Hall Temp Difference (Tudiff) 
СО ет? 


! CALCULATE LOCAL CHANNEL TEMPERATURES AT INLET. MIDDLE AND OUTLET. 


Toin=Tint(.S*Tdiff/12) 
Tomid=Tint(4.5+*Tdiff/12) 
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= 30. 


Тсоит= ТПО БТЕ ИАС 

PRINT USING "12X,""Average Local Channel Temp at Inlet (Tcin) = єр. 
СОО "ОТЕТ 

PRINT USING "12X,""Average Local Channel Temp at Middle (Tcmid) = “р. 
СС" " етіс 


PRINT USING "12X,""Average Local Channel Temp at Outlet (Tcout) = "ИЗ. 
(C)*""2mmout 


! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT INLET. 


Sum=0. 
Iadd-0 
IF Ipos-1 THEN Iadd-40 


! CORRECTION FOR BROKEN T/C CHANNELS 5, AND 47 CURVED HEATER ONLY. 


IF Itype-1 THEN 

РОК 1=0 Т0 6 

IF (1Іроѕ=0 АМ) І-5) OR (Ipos-1 AND I-6) THEN 3810 
Sum*Sum*T(I*Iadd) 

NEXT I 

Twin=Sum/6 

Se 


Е 

FOR I=0 TO 7 

IF (Ipos-0 AND I-7) OR (Ipos-! AND I=6) THEN 3870 
Sum*Sum*T(I*Iadd) 

NEXT I 

Тшіп"бит// 

END IF 


PRINT USING "12X,""Avg Local Heated Hall Temp at Inlet (Tuin) ===" 3D. 
CC ete ІШІПТІ 


! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE NEAR MIDDLE. 
Sum=0. 


! CORRECTION FOR BROKEN T/C CHANNEL 16 CURVED HEATER ONLY. 


IF Itype=1 AND Ipos=0 THEN 
РОК І=9 ТО 15 

ТЕ 1=10 ТНЕМ 4020 

Sum=Sum+T (1) 

NEXT I 

Twmid=Sum/6 

ENSE 

FOR 1=9 ТО 16 

IF (Ipos=0 AND I-10) OR (Ipos=1 AND I-9) THEN 4080 
Sum*Sum*T(I*Iadd) 

NEXT I 

Twmid=Sum/7 

END IF 

PRINT USING "12X,""Ava Local Heated Hall Temp Near Middle (Tumid) = "",3D. 
Пејн Тота 


| CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT OUTLET. 
Sum=0. 


! CORRECTION FOR BROKEN T/C CHANNEL 26 CURVED HEATER ONLY. 
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IF Itype=1 AND Ipos=0 THEN 
FOR 1923 ТО 29 

IF I-26 THEN 4230 
Sum*Sum*T(I) 

NEXT I 

Twout*Sum/6 

CUSE 

FOR I=23 TO 29 
Sum*Sum*T(I*Iadd?) 

NEXT I 

Tuout*Sum/7 

END IF 

PRINT USING "12X,""Avg Local Heated Hall Temp at Outlet (Tuout) „р. 
(С)”“” Ковић 
PRINT 

PRINT 


! ENTER ORIFICE PRESSURE DATA. 


BEEP 


INPUT "ENTER PATM (inHg), DPM(inH20). P1M(inH20), REI°.Patm.Dpm,P1m.Rel 
PRINT USING °SX,"" THE FOLLOWING ORIFICE DATA WHERE ENTERED:""" 

PRINT 

PRINT USING "12X,""Patm(i1nHa) DPM(1nH20) P1M(1nH20) REI""" 
PRINT USING °14X,2(2D.0D,10X),(2D.DD,7X) ,S5D.D" ;Patm,Dpoem,Pim.Rel 

PRINT 


! CONVERT PRESSURE READINGS TO SI UNITS AND P1M TO ABSOLUTE. 


Patm=Patm#33/76.8 
Pdel=Dpm*248 .84 

Р1 =Ра+т- (Р1т*248.84) 
ТогҒ“Т(33) 


! PHYSICAL PROPERTIES AND CONSTANTS. 


R=286.987 !FOR AIR. 
Rho=sP1/(R*(Torf+273.15)) !DENSITY OF AIR AT ORIF PLATE. 
Gamma=1.40 !FOR AIR. 


ЕР ІСр GOOD FOR Tblk BETWEEN 12 AND 33 (DEG C). 
с=1. 


51ата=5 .6698Е-8 
Еси«.12 ‘EMISSIVITY OF COPPER. 


! DIMENSIONS OF CHANNEL. 


De=.006350 
Wid=.2540 
Ri=.297 

Рше+=2* (0с+Н19) 
Ас=0с*Н 19 
Dhd=4*Ac/Pwet 
Apl=.0774192 


PROPERTIES OF INSURATION. 


X1ins*.0127 
Kins*3.8E-2 


! ENTER ORIFICE CONFIGURATION. 


87 


4780! 
4790 
4800 
4810 
4820 
4830! 
4840 
4850! 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
49401 
4950 
4960! 
4970 
4980! 
4990 
5000! 
5010 


ВЕЕР 

INPUT “SELECT DIAMETER OF ORIFICE (0=.5334, 1t=1.0755 (inches))",Size 
IF Size=0 THEN Dorf=.013548 

IF Size=1 THEN Dorf=.027318 


! DIMENSIONS OF PIPE AND ORIFICE PLATE. 
A-(PI*Dorf'2)/4 


Dpipe=.051772 
Apipe-PI*Dpipe'2/4 


Beta=Dorf/Dpipe 

PRINT USING "12X,""Dorf (m) A(m*2) Beta""* 
PRINT USING °12X,(Z.6D,6X) ,(Z.3DE,6X) ,Z.4D";Dorf,A,Beta 
PRINT 

PRINT 


! CORRELATION FOR EXPANSION FACTOR BASED ON 1D AND 1/2D TAPS. 
Y=1-(.41+.35*Beta’4)*(Pdel/(Gamma*P1)) 

! CALCULATION OF FLOW COEFFICIENT (K) IN SI UNITS. 

B= .0002+2.794E-S/Dpipe+( .0038+1.016E-S/Dpipe)*(Beta*2+(16.5+1.968504E+2*Dp 


ipe)*Beta' 16) 


5020! 
5030 
5040! 
5050 
5060 
*Beta' 
5070 
5080 
5090! 
5100 
5110! 
5120 
5130! 
5140 
5150! 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 


304%: 


201 
5270! 
5280 
5290! 
5300 
5310 
5320 
5330! 


! КТ АМ К2 АКЕ DUMMY VARIABLES TO CALCULATE Ko. 


K1-.6014-5.3974064E-3/Dpipe'.25 
K2-(.376*2.8971138bE-2/Dpipe'.25)*(1.6129E-7/(Dpipe'2*Beta'2*6.35E-5*Dpipe) 
4*1.5*Beta' 16) 

Ко=К1 +К2 

К=Ко+1000*В/Ке1 ^.5 


! Mu CORRELATION GOOD FOR Torf BETWEEN 17 AND 41 DEGREES CELSIUS. 
Mu=4.6971429E-8«Torf+1.7194722E-5 
! CALCULATE Mdot AND REpipe AND COMPARE TO PREDICTED КЕ1. 


Mdot-K*A*Y*x(2*Gc*Rho*xPdel)'.5 

Repipe=(Mdot*Dpipe) /(Mu*Apipe) 

Dif f=(ABS(Rel-Repipe) /Rel)*«100 

Rel=Repipe 

IF Diff>.001 THEN 5080 

EE USING “9X,"“"“THE FOLLOWING DATA HERE CALCULATED: *"" 

PRINT USING "12X,""Orifice Expansion Factor (Y) = 44. 

PRINT USING "12X,""Orifice Flov Coefficient (К) - і 

г USING "12X,""“Density Based on Torf (Rho) = 744% 
о 

PRINT USING "12X,""Viscosity Based on Torf (Mu) = 73027 ЗВЕ. "“ (сайп е) 


(ка/п” 


! Mu CORRELATION GOOD FOR Tblk BETWEEN 17 AND 41 DEGREES CELSIUS. 
Mu-4.6971429E-8*Tb1k *1.7194722E -5 


Red=(Mdot*Dc)/(Nu*Ac) 
Rehd=(Mdot*Dhd) / (Mu*Ac) 
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5340 
5350! 
5360 
5370 
5380 
5390 
5400 
5410! 
5420 
5430! 
5440 
5450! 
5460 
54701 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570! 
5580 
5590! 
5500 
5610! 
5620 
5630! 
5640 
5650 
5660! 
5670 
56801 
5690 
5700 
5710 
5720 
5730 
5740 
5750 
5760 
5770 
5780! 
5790 
5800! 
5810 
5820 
ни Но 
5830 
ека 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 


! CALCULATION OF POMER INTO HEATER PLATE. 


Tadd=0 

IF Ipos=1 THEN Iadd=39 

IF Ipos-0 THEN Rpr=2.01729 

IF Ipos=1 THEN Врг«2.00839 

Qp=(Emf (38+1ladd) *Emf (39+Iadd))/Rpr 


! CALCULATION OF HEAT INPUT TO AIR. 
Qair=Mdot*Cp*( Tout- Tin) 
! CALCULATE LOCAL HEAT TRANSFER COEF. AND NUSSELT NUMBER. 


Kairz7.7257143E-5*Tcin*.024165836 
Н. п=0аіт/ (Ар 1% (Тшіп-Тсіп) ) 
Nuin=(Hin*Dhd) /Kair 
Kair-7.7257143E-5*Tcmid*.024165836 
Hmids*Qair/(fipl*(Tumid-Tcmid)) 
Numid=(Hmid*Dhd) /Kair 

Каіг=7. 7257 143Е-5 + Тсоиї+. 024165836 
Нои = 0а іт / (Ар 1 * (Тшоиї- Тсоиї ) ) 
Nuout*(Hout*Dhd) /Kair 


! Kair CORRELATION GOOD FOR Tblk BETWEEN 17 AND 41 DEGREES CELSIUS. 
Kair=7.7257143E-S*Tblk+.024165836 
! CALCULATE AVERAGE HEAT TRANSFER COEF. AND NUSSELT NUMBER. 


Намаєда1 г / (Ар 1 ж Тае) 
Nuavga* (Hava*Dhd) /Kair 


! CALCULATE HEAT LOSSES. 


Qlol-*(Kins*fpl*(Tinsol-Tinso2))/Xins 
Qlo2-(Kins*fipl*(Tinso2-Tinso3))/Xins 

010= (0101+0102) /2 
Qlil=(Kins*Ap1*(Tinsil-Tinsi2))/Xins 
Q1i2=(Kins*Apl*(Tinsi2-Tinsi3))/Xins 
011=(0111+0112)/2 
Кг“2.0%((1-Еси) / (ЙріжЕси))%(1.0/Ар1) 

Qr=(Sigma* (ABS (( Twot273.15) *4-(Twit+273.15)°4)))/Rr 
Qdel-Qp-Qair-Qlo!-Qr-Qli! 


! CALCULATE DEAN NUMBER. 


De=.S*Red*((Dc/2)/Ri)*.5 
PRINT USING “12X,""Viscosity Based on Tblk (Mu) = "7. 22.308.“ кан 5) 


PRINT USING "12X,""Therm Cond Based on Tblk (Kair) = ^"",Z.3DE."" (M/m.K)* 
r 

PRINT 

PRINT USING "12X,""Patm x "".6D.DD."" (N/m'2)""":Patm 
PRINT USING "12X,""Pdel = + роб == сМутог)“"" Рвеј 
PRINT USING "12X,""P1 ОВО СМУ те Е 
PRINT USING "12X,""Mdot = и. (kg/s) S dot 
PRINT USING "12X,""Repipe = "",6D.D":Repipe 

PRINT USING "12X,""Red = "" 6D.D":Red 

PRINT USING "12X,""Rehd » "".6D.D":Rehd 
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5920! 
5930 
5940! 
5950 
5960 
5870 


5980 


! ONLY PRINT DEAN NUMBER WHEN HEATING CURVED SECTION. 


IF Itype-1 THEN PRINT USING "12X,""De = "*,4D.D";De 
PRINT 


PRINT USING "12Х," "др « "Є 30,30,"" (Watts) "" 7X, 241215 202212 5,309 (A 
atts)"""<@e7Ga7r 


PRINT USING *12Х," "0101 = "",2D.3D,"* (Matts)"*,7X, ""Qli1 


(Нас іє) 77979 До | 
5990 PRINT USING "12X,""01o2 <. "*,2D.3D,"'" (Matts)"",7X.""Q01i2 
(Natts)""*":0102,Q11i2 


6000 


PRINT USING "12X,""Qlo 


"*.,2D.3D.,"" (Matte) "SX НАШЕ 


(Наі%5)"”"%0 10 ВІ 


6010 PRINT USING "12Х, ""Ог = ”"”(20:3 2 “2” (На 7”"”,7Х Обе! 

(Watts)""";Qr,Qdel 

6020 PRINT 

ш PRINT USING "12Х,""Ніп = "*,3D.DD,"" (H/m'2C)"*",7X, ""Nuin 

in,Nuin 

6040 PRINT USING "12X,""Hmid < "",3D.DD,"" (W/m°2C)"",7X,""Numid 
mid,Numid 

6050 PRINT USING "12X,""Hout » "",3D.DD,"" (W/m'2C)"",7X,""Nuout 


out,Nuout 


6060 


PRINT USING "12X,"°Havg = "",S3D.DD,"" (W/m'2C)"* "*,7X, " "Nuava 


avg,Nuava 


6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210! 
6220 
6230! 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320! 
6330 
6340! 
6350 
6360 
6370 
6380 
6390 
6400 
5410 
6420 


PRINT 

PRINT 

BEEP 

INPUT “WILL THERE BE ANOTHER RUN? (1<Ү,0- №)“, Со оп 
IF бо оп<>0 ТНЕМ 

PRINT USING "&,»" 

IF Im=0 THEN WAIT 600 

GOTO 1700 

END IF 

PRINT USING *8X."*END OF RUN"*"*" 
PRINT USING *&,*z" 

PRINTER IS 1 

ASSIGN @File TO * 

END 


! CONVERTS EMF TO DEGREES CELSIUS. 


SUB Tvsv(V,T) 

COM /Со/ 0(7),Ва(76,2) 
Sum=0 

FOR I-0 TO ? 
Sum=Sum+D( I) *V°TI 

NEXT I 

T=Sum 

SUBEND 


! CALIBRATES T/C READINGS. 


DEF ЕМТсогг(Т,Г) 

СОМ /Со/ р(7),Аа(76,2) 
Тс=Аа (1,0) 

FOR J=1 TO 2 

Тес“ Тс+ћа(1,Ј)+Т Ј 

МЕХТ Ј 

RETURN Tc 

FNEND 
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APPENDIX B: SAMPLE CALCULATIONS 


Figure 36 below, shows the major heat transfer 
components for each of the test sections. The sample 
calculations that follow, demonstrate the methods used by 
the computer to calculate these components as well as the 
Reynolds number, average heat transfer coefficient, and 
average Nusselt number for each set of data. Part A of 
Appendix B contains a sample printout of a data run, 
followed by a summation of the data. The sample 
calculations are for the concave curved section, but those 


for the straight section are similar. 


HEATED COPPER WALL 


Qlo 





Figure 36. Energy Balance in Test Section. 
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1. SAMPLE CALCULATION DATA 


THE FOLLOWING CHANNEL CONFIGURATION WAS SELECTED: 


Heating CURVED Test Section. 
Heating OUTER Plate. 


Month, Day, and Time: 


THE FOLLOWING DATA WERE RECORDED: 


Heater Voltaae (Vh) 
Resistor Voltage (Vr) 


Precision Resistor (Rpr) = 


OUTER PLATE; 


T/C Number: 


Temp (С): 


Т/С Number: 


Temp (C): 
T/C Number: 
Temp (C): 
T/C Number: 
Тетр (С): 
T/C Number: 
Temp (C); 


OUTER INSULATION: 


T/C Number: 
Тетр (С): 


ORIFICE TEMP (Torf) = 24.16 (С) 
INLET TEMPERATURE: 
T/C Number: 35 36 
Temp (C): 202255 20.48 
INNER PLATE: 
T/C Number: 41 42 
Temp (C): 21.43 217251 
T/C Number: 47 48 
Тетр (С): ОЗ 22753 
[/C Number: 53 54 
Temp (C): оно 21.27 
Т/С Number: E] 60 
Temp (C): DS 217556 
T/C Number: 65 66 
Temp (C): 22.05 212 
INNER INSULATION: 
T/C Number: 71 (Tinsi 1) 
Temp (C); 21253 
QUTLET TEMPERATURE: 
T/C Number: 74 ЙЕ 
Temp (C): 24.29 24.38 


25:605 337 20 
= ЭБ CU 


46.47 


31 (Tinsol) 
е: 


92 


l2 


US 


ПО ЗА 


2.01729 (Оћт5) 


4 

дно 
10 

43.47 
16 


41.97 
22 

33:65 
28 

45.80 


32 (Tinso2) 


37.29 
37 38 
20.47 20.43 
43 44 
21.16 gr 
49 50 
20. Dn 21.05 
55 56 
мерни 21.19 
61 62 
21.47 21.74 
67 68 
21.85 ПИВ 
(Tinsi2) 
211221 
76 77 
24.38 24.25 


3 
44,55 
11 
44.77 
Í 


1 
38.06 

Ze 
45-32 

23 
дого 


(Tinso3) 


23835 


45 
21715 
5 | 
ИЗО 
-2 
2115 
63 
21855 
69 
21.74 


(Tinsi3) 


еу; 


46 
21212 

52 
21.46 
21.35 

64 
2211 
2113 


THE FOLLOWING TEMPERATURES WERE CALCULATED: 


Average Outer Wall Temperature (Two) 
Average Inner Wall Temperature (Twi) 
Average Outlet Temperature (Tout) 
Average Inlet Temperature (Tin) 


Channel Inlet and Outlet Temp Diff (Tdiff) 


Average Bulk Temperature (Tblk) 
Mean Temperature Difference (Tdel) 


Outer and Inner Wall Temp Difference (Tudiff) 
Average Local Channel Temp at Inlet (Tcin) 
Average Local Channel Temp at Middle (Tcmid) 
Average Local Channel Temp at Outlet (Tcout) 
Ava Local Heated Hall Temp at Inlet (Tuin) 

Avg Local Heated Wall Temp Near Middle (Twmid) 
Avg Local Heated Wall Temp at Outlet (Tuout) 


THE FOLLOWING ORIFICE DATA WERE ENTERED: 


Patm( inHg) DPM(1nH20) 


29.79 18.58 
Dorf(m) A(m*2) 
0.027218 5.861Е-04 


THE FOLLOWING DATA WERE CALCULATED: 


Orifice Expansion Factor (Y) 
Orifice Flow Coefficient (K) 
Density Based on Torf (Rho) 

Viscosity Based on llorf (Mu) 
Viscosity Based on Tblk (Ми) 


Therm Cond Based on Tblk (Kair?) 


Patm - 100594.87 (N/m'2) 
Pdel = 4623.45 (М/т 2) 
P1 - 974/71.93 (М/т 2) 
Mdot - 0.0375 (kg/s) 
Кер:ре = 50352.2 

Red = 8097.1 

Rehd .-.. MEI 

De = 418.6 

Ор “ПЕ, ӘЛІ Watts) 
0101 = 2.847 (Watts) 
Qlo2 = 1.838 (Watts) 

Qlo = 2.342 (Watts) 

От = СІ Ма ска) 

Hin = 81.90 (W/m*2C) 
Hmid = 86.78 (W/m*2C) 
Hout = 85.51 (И/т* 2С) 
Нама =  83.75 (М/т 20) 


9,3 


= 44, 
x 2], 
. 24, 
а, 20 
= E 
= 
чо 
= 23 
= 20 
= 21 
= 24 
= 43 
= 43. 
= 46. 
Р1МСїпН20) RE! 
12.55 15500.0 
Beta 
0.5277 

= 0.9852 

= 0.6323 

= 1.1424 (Ка/т 3) 

- 18.330E-06 (kg/m.s) 

- 18.247E-06 (kg/m.s) 

= 2.590E-02 (W/m.K) 


Ја 
0111 
0112 
011 

Оде 1 


Nuin 

Numid 
Nuout 
Nuavg 


noe’ uw ott 


144.950 (Watts) 
ООЗЕ B3 tts ) 
.w038 (Hatts) 
.023 (Watts) 

9.415 (Watts) 


33725 
41.58 
40.70 
40.07 


ge 
(C) 
(С) 
або 
(С) 
(C) 
(C) 
(С) 
(C) 
(C) 
(C) 
(C) 
(C) 
(C) 


V = 61.337 У 


Н 
Vpn = 5. 195 Ч 
Кър = 2.0173 RQ 
p = 29.798 n. На 
AP = 18.58 in. Но 
Р, ш 12,554 п. H,O 
ре = 20.48 С 
in 
TOC = 24.32 С 
- = 44.76 C 
: = 21.49 С 
wi 
DIEN ЕП . 53 16 
Tins,i2 = 21.4908 
Ti nemo = 49.58 С 
Tins,02 = Ey 29 C 
Torf =m, 16 
_ -3 
Ko 2258009 10 W/mC 
Сраїг = 1006 Јука К 
цэс z 4.6971 x 10 X Tm + 1.71947 x 10 
K. 3.80 x 10 ^ W/mC 
ins 
AXi ns 20-012 78m 
€ = 0.12 
си 
0 NES 10-17: 
В =e 52777 
Y = 1.40 
2 
gc - ] Kg m/Ns 


94 


мо-мі 


286.98 Nm/Kg K 


= 1.0 


. 07742 т^ 


. 0016 m 


„00211 п? 


‚00635 т 


‚0518 т 


‚027318 м 


25207m 


. 297 т 


2. TEMPERATURE CALCULATIONS 


a. 


Bulk Temperature (T lk? 


T . in (| out _ 28.48 % 21:22 = 22.40 С 


Mean Temperature Difference (AT) 


АТ = Т M = 44.76 - 22.40 = 22.36 


wo blk 


Temperature Difference Between Outer and Inner Wall 


шэг. 


т. 22214276 /-2721 09"75:423.27-С 
wi 


КЕСТІ Too 


95 


3. | POWER CALCULATIONS 


а. 


Power Supplied (0Q.) 


P 


-R gels si oS) № 
Ор ES 12-0173] - 157.97 W 


Heat Lost Through the Back of the Outer Plate (Qi)? 


lo -2 
(АХ, 5) / (Кур) (Арг) (02-0127)/ (3.80х10 _) 6. 07 о 
= 2.85 М 


Heat Lost Through the Back of the Inner Plate (01) 


E SEE MIS 2 (21.53 2 р 


Q E нн а D ааа 
1 -2 
бах, <) И (К; 2) (Арр) (0.0127)/(3.80х10 °) (.07742) 


1 


2010 


Heat Radiated (Q,) 


(1) Radiation Resistance (Rp) 





1 - *cu 1 
В = 2х (—) + 
R Apr * cu Apr Е 20-11 
1 2 1 2 -2 
-- г - 1) = ————— [—- 1] » 202.36 m 
Арт E eH 0.07742 2152 


(2) Heat Radiated (Q,) 


4 4 


c (T - CT. ) 
о, E No wi T [=] K 
R 
-8 4 à 
- (5.669x10 7)[(44.76*273)" - (21.49*273)"] 
2 02. 36 
= 0.75 И 


4. MASS FLOW RATE CALCULATIONS 
a. Pressure Conversions 


(КГР 29.79 in Hg x 3376.8 - 100594.87 N/m* 


atm 
(2) АР = 18.58 in Н.О х 248.84 = 4623.5 N/m* 
(3) P, = (12.55 іп Н.О х 248.84) - 100594.87 


97471.93 м/п? 


tl 


air 
EET (97471.93) 3 
2-1: т койш (286.98) (24.16 + 273) АИ 
c. Expansion Factor (Y) 


EUM cs. 
ҮР, 


4 4623.5 
1 - [.41 + .35(.5277) ] 1.40) (97471.93) 


. 9852 


B 


d. Area of Orifice (A) 
0.005 о : 
А = —————— = ———— 2,005861 m 
4 4 

e. Mass Flow Rate (m) 

m= YKA 2g, o.,, AP = (.9852)K(.0005861) 

x 2 (1) (11424) (462375) 

m = .05935 K 

Iterating: 


Assume a Pipe Reynolds number = 50000 
Obtain a value for K, the flow coefficient, from 


reference 30. К = „6325 


Solve for m. m = .0375 Kg/s 


- P Upipe . — (.0375)(.0518) | 


A (400210 ) (1893351088 


Solve for new Re pipe 
рїре"аїү 


Б - 550225 


е_. 

ртре 
Check convergence and repeat process if necessary. 
(Convergence if difference less than .001) 


m = .0375 Kg/s 


98 


REYNOLDS NUMBER CALCULATIONS 


m 
ыг RED i pe = pipe _ (:03757 000518) E 50225 
dt: 0800211) 01823320009) 
pipe air 
m D 
ж с _ _(.0375) (.00635) _ _ sog; | 
Aa Magy — (40016) (18. 25x10 9) 
А оп Dag 
Бі» *hd A yu. 
© alr 
4x A 
С (4) (. 0016) 
ГЕ | О 01229 п 
hd Poo. 75207 
m Dia (110221 
Вера = нн = 15799.1 
UL ION OUT 25,510) 
HEAT CONVECTED TO AIR CALCULATION 
а. 9,1,2 ПС,(Т,,-Ті,) 5 (20375) (1006) (24.32-20.48) 2144.95 W 


AVERAGE HEAT TRANSFER COEFFICIENT CALCULATION 


Галат 134-95 EN - 55. 72 AC 


99 


9. DEAN NUMBER CALCULATION 


= 418.6 





Rea  /pc/2 _ 8097.1 .00635/2 
2 ПИНО и 


даж Ri „297 


100 


APPENDIX C:  EXPERIMENTAL UÜNCERTATNTY 


The uncertainties for the major variable in the 
experiments, were calculated in accordance with the method 
described by S. Kline and F. McClintoch [Ref. 33]. The 
estimates of the uncertainty in the measured quantities were 
made conservatively. А5 а result, there is considerable 
confidence in the uncertainties as calculated. 

The following equations were used to calculate the 


uncertainties: 




















1/2 
2 ар, 2 ат 2 
ар _ ав 1 ої 
(1) Бағы (Ro + + (F ) ) 
I orf 
: 2 2 2 
ам ах ак дА 
(2) Шо КИ хс) КЕ) + aa ) 
m 
272 
ар: 2 2 
alr dAP 
* (1/4 Mc t Gp)? ) ] 
air 
1/2 
аф іг Ж dm 2 а Coair 2 ИЕС 102 с 
(221) сл л Мо + бус) + т ТТ ) ) 
а1г m pair out in 
172 
-> ао: аА 2 2 
аһ аіг PL AT 
(4) => = (( Ји et } м) ) ) 
h Qair А рі, A 
m 1/2 
dNu = ар 2 dK.. 2 
(5) ——%=- (197+ (,99) + (=) ) 
Мына ћ ћа alr 


ТОД 


dRe, 4 Р 


Вега 


ар 


(6) 5 





52222 
(олг 
m 


The following quantities 
calculated by dividing the es 
by the value of the quantity, 
considered constant for the r 


experiment: 


Quantity 


A 

Ас 
АРІ, 
Сраїг 
рс 


ћа 


а1г 


atm 
1а 
alr 


The following quantities 


hydraulic Reynolds numbers be 





hd 


2 


2 
па 
) 


ОН 
) + 


T 





алар 


had their uncertainty 
timated error in the quantity 
and these uncertainties are 


ange of values in this 


Uncertainty 
. 0093 


. 0408 
. 0080 
. 0020 
. 0400 
320/596 
.0030 
. 0038 
‚0010 
„0003 
. 0080 
. 0020 


. 0054 


had the worst case error for 


low 17000: 


102 


Quantity Erron 


Tblk в.17с 
Тёре зале 
orf %.17с 
ЛЕК iG 
Dus +.1°с 
Tout lin %.27с 
АТ +. 226 


For hydraulic Reynolds numbers greater than 17000 the 
error in the outlet temperature was much higher, as 


mentioned previously, and the following errors apply: 


Quantity Error 
Tolk +.2°C 
T С 
I РА С 
Tout %.37с 
Tos а 
Tout lin %.47с 
АТ +.3°с 


The uncertainties for the following quantities were 
calculated using the sample data for the outer curved test 


section at a hydraulic Reynolds number of 15800: 


Коз 


Quantity Uncertainty 


Р, - 0159 
АР .0108 
Torf . 0003 
Tout en • 0521 
АТ . 0089 
aee . 0159 
m „0139 
Ое . 0540 
h 70953 
Мира . 0681 
Кера . 0588 


The uncertainty for the following quantities were 
calculated for the outer curved test section at a hydraulic 


Reynolds number of 21200: 


Quantity Uncertainty 
Р, „0043 

АР . 0053 
TS .0003 

AT „ОЛА 
ЕРТЕ .0043 

m . 0105 
Deu .1044 


104 


h .1053 
Ми #7126 


“0581 


By multiplying the quantity by its uncertainty the 
possible error can be calculated. The previous uncertainty 
calculations for a hydraulic Reynolds number of 15800 
yielded Nu,4, - 40.07 *2.73 and a Re,4 - 15800 + 928. For 
the run at a hydraulic Reynolds number of 21200 the Ми, з = 
44.71 + 5.03 апа а Ке а ЕВРО 1232. 

The major source of uncertainty in the average Nusselt 
number is the uncertainty in the difference between the 
outlet and inlet air temperature. If this uncertainty can 
be brought below four percent, then the uncertainty of the 
average Nusselt number will be dictated by the uncertainty 
of the measured dimensions of the channel. In particular, 
the channel height (Dc) which has an uncertainty dD./D. = 
.0400. The uncertainty of the channel height can only be 


reduced by increasing the channel height, and/or using 


materials with better tolerances. 


105 


12972 


LISTESEOREBERENGES 


Taylor, б. I., "Stability of a Уїзссйз 2 Ё (1! Ээ 
Between Two Rotating Cylinders," Philosophical Tran- 


sactions of the Royal Society of London, series A, 
V. 233, рота. 


National Advisory Committee for Aeronautics, Technical 


Memorandum 1375, On the Three Dimensional Instability 


of Laminar Boundary Layers on Concave Walls, by H. 
Gortler, 1942. 


Smith, A. M. O., "On the Growth of Taylor-Gortler 
Vortices Along Highly Concave Walls," Quarterly of 


Applied Mathematics, V. 8, рр. 233-263, November 1955. 


Dement'eva, K. V. and Aronov, I. Z., "Hydrodynamics and 
Heat Transfer in Curvelinear Channels of Rectangular 


Cross Section," Journal of Engineering Physics, V. 34, 
No. 6, pp. 666-671, 1978. 


Mayle, R. E., Kopper, F. C., Blair, M. F., and Bailey, 
D. A., "Effect of Streamline Curvature on Film 


Cooling," Journal of Engineering for Power, Trans. 
ASME, V. 99, Series A, No. l, pp. 77-82, January 1977. 


Nicolas, J. and LeMeur, A., "Curvature Effects on a 


Turbine Blade Cooling Film," ASME Paper No. 74-GT-156, 
1974. 


Folayan, C. O. and Whitelaw, J. H., "The Effectiveness 
of Two-Dimensional Film-Cooling Over Curved Surfaces," 


ДЭМЕ-Рарег Мо. 20-НТ 310 Ge 


Lord Rayleigh, "On the Dynamics of Revolving Fluids,' 
Proceedings of The Royal Society of London, Series A; 
V. 93, pp. 148-154, 1916. Reprints in Scientific 
Papers, V. 6, pp. 447-453. 


Taylor, G. I., “Distribution of мена rand 
Temperature Between Concentric Rotating Cylinders," 


Proceedings of the Royal Society ошшканосн, Series A, 
У. 151, pp. 494-512, 9 


Dean, W. R., "Fluid Motion in a Curved Channel," 


Proceedings of the Royal Society of London, Series A, 
V. 121, pp. 402-420, 1929: 


106 


ШЕЛ 


ї 2% 


ЫЗ, 


14. 


ts: 


1862 


т 


I8. 


INO 


20. 


212 


Reidi "W. H., “On the Stabildty- of VisScous Ећом- 1п га 
Curved Channel," Proceedings of the Royal Society of 
London, Series A, V. 244, pp. 186-198, 1958. 
Schlichting: H Boundary т Theory, 7th Ed., 

pp. 529-536, McGraw-Hill, 1979. 

Kelleher, M. D., Flentie, D. L., and McKee, R. J., "An 
Experimental Study of the Secondary Flow in a Curved 


Rectangular Channel," Journal of Fluids Engineering, 
М2 02. "реке 92 96March 1980. 


ОО Т 20017000. F. G., and Crane, R. I., 
"Measurement within Gortler Vortices," Journal of 
ЕТО ШЕП С СЛОИ Е ОШ pp. 517-520, December 1979. 


Aihara, Y., "Nonlinear Analysis of Gortler Vortices,"l 
The Physics of Fluids, 19, pp. l655-1660, November 
1976. 


Kreith, F., "The Influence of Curvature on Heat 
Transfer to Incompressible Fluids," Trans. ASME, V. 77, 
рр. 1247-1256, 1955. 


AemospaceMkesearch Laboratories Reprt ARL 65-68, А 
Simplified Approachwto the Influence of Gortler-Type 
Vortices on the Heat-Transfer from a Wall, by Leif N. 
Persen, May 1965. 


Cheng, K. C., and Akiyama, M., “Laminar Forced Convec- 
tion Heat Transfer in Curved Rectangular Channels," 
International Journal of Heat and Mass Transfer, V. 13, 
pe. 471-490, 1970. 


Shibani, A. A., and Ozisik, M. N., “A Solution to Heat 
Transfer in Turbulent Flow Between Parallel Plates," 
International Journal of Heat and Mass Transfer, V. 20, 
БЕН 65-072, 1977. 


Mori, Y., Uchida, Y., and Ukon, T., "Forced Convective 
Heat Transfer in a Curved Channel with a Square Cross 
Section," International Journal of Heat and Mass 
Hamster. V. 14, рр. 1787-1805, 1971. 


Kays, W. M. and Leung, E. Y., "Heat Transfer in Annular 
Passages - Hydrodynamically Developed Turbulent Flow 


with Arbitrarily Prescribed Heat Flux," International 
Journal of Heat and Mass Transfer, V. 6, pp. 507-557, 


ОР 


107 


22: 


23: 


24. 


258 


26% 


275 


28. 


298 


30i 


31 


а 


33% 


Brinich, Pæ F. and Graham Rə W., "Flowwandaldgest 
Transfer in a Curved Channel," NASA Technical Note No. 
TN-D-8464, 1977. 


McKee, R. J., An Experimental Study of Taylor-Gortler 


Vortices in a Rectangular Channel, Eng. Thesis, Naval 
Postgraduate School, Monterey, California, 1973. 


Durao, M. do Carmo, Investigation of Heat Transfer in 
Straight and Curved Rectangular Ducts Using Liquid 


Crystal Thermography, Eng. Thesis, Naval Postgraduate 
School, Monterey, California, 1977. 


Ballard, J. C. III, Investigationwoft Heat Transfer in 


Straight and Curved Rectangular Ducts, Master's Thesis, 
Naval Postgraduate School, Monterey, California, 1980. 


Holihan, R. G., Jr., Investigation of Heat Transfer in 
Straight and Curved Rectangular Ducts for Laminar and 


Transition Flows, Master's Thesis, Naval Postgraduate 
School, Monterey, California, 1981. 


Daughety, S. F., Experimental Investigation of 
Turbulent Heat Transfer in Straight and Curve 
Rectangular Ducts, Master's Thesis, Naval Postgraduate 
School, Monterey, California, 1983. 


Wilson, J. L., Experimental Investigation of Turbulent 
Heat Transfer in Straight and Curved Rectanglar Ducts, 
Monterey, 


Master's Thesis, Naval Postgraduate School, 
California, 1984. 


Gebhart, B., Heat Transfer, 2nd Edition, p. 260, 
McGraw-Hill, New York, 1971. 


The American Society of Mechanical Engineers, 
Supplement to ASME Power Test Codes, Chapter 4, Flow 
Measurement, p. 25, 1959. 


Benedict, R. P., Fundamentals of Temperature, Pressure, 
and Flow Measurements, 2nd Edition, pp. 443-447, John 
Wiley and Sons, 1969. 


Shah, R. K. and London, A. L., Laminar Flow Forced 
Convection in Ducts, Supplement 1, pp. 305-312, 


Academic Press, 1978. 


Kline, S. J., and McClintock, F. A., "Describing 
Uncertainties in Single-Sample Experiments," Mechanical 


Engineering, V. 75, pp. 3-8, January 1953. 


108 


ШТАТ DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22304-6145 


Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943-5100 


Department Chairman, Code 69 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5100 


Professor M. D. Kelleher, Code 69Kk 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93943-5100 


Lieutenant George G. Galyo 
724 Charette Road 


Philadelphia, 


Pennsylvania 19115 


109 


No. Copies 
2 











Thesis 
С1422 


GL 





216455 


Galyo 
Experimental inves- 
tigation of turbulent 
heat transfer in 
straight and curved 
rectangular ducts. 











~ LN m EN E 
> ~ - = е - ~ ЕЗ - ~ + ~ ч. 
- - > ГЭ am BOS = - ~ ть 
E -- “ы қ - сэ - ere 5 Ld „= ©, “ - » з қ, 
- – - ~ - .... = ~ БЭ -- ГЭ ` 
т ..... - - ~ LS е ~ ~ ч ~ . + 
.... - ~ P P ан ~ ~ ove У эт. - - Ы « Ы <=. . А М най > 
E ре ы ~ - РЕЗЕ & ~ ~ БД - E - — 2 úm - ч ме ын = » 
E —! жеты ~ - ~ ~ „= <» == = E “ - хү he > "о = - - у, 9 eec 5 - Б 
м^ -.. 2 јр PE E - = * ы ~ ~ ~ ГИ - 2 E La ~ e P ға” а ~ ~ = ле = pu a = 
n ~ ~ = == - ++ РИК = ы ~ ~ = ы E ы ж - ara s е - ” = 2 up ЭЭ 29 s a 
— - Р ~ -- ~ ~ E ^n - “ 3 P 4. т ы ө ~ + = + - - = e + “и — - P - та? a E 
E = — M ње и ы ~ ~, - аль м я E - - - я E P e e o - ——- => sma. — 2 -. 
= © ъ= A „= -- - d» e oc “> - ~ жы % E ы ы 7 тез” РК = - "P E 4 T 
- rt te] T п—— - 0t ——Ó—— = - шог 4. еа ? A i» е ы = ы - s E = = ы “= pp par» = -"” .” sad adiens n Pe 
Mad Xs lind a PM m ee чеч E d - "X on wv ~ æ ~ е ° = "У "m - P - ^. "DP - мэ”? = P оо ль о А е о ә 
-— . > ~ ~ pu > ч ч ~ а "TL d Б ~ = "E = o E Ын = тет ә =» - = a Ф = - д“?, рл” у" oe Е ПР чи ы m = 
Р - `% ~- й E ” -— ^ —^ Б " ns * e ~ чы, sar ТТ 2? ~ те Ы > ~ f Leo я е 2 Ч - E ^ = ш ы Б ~ 2 =“ ° = E == e - 
мэ - - om cmd Т е 12 - ~ - ы ^ oo - >: - ын = €— = ыы ИЙ ва сх Led t ий. > ~ тээж» - - oon - а «әс 
be P РСЕ ч » — n ... ~ PN ч E ~ • г I ^ - > ~ . Р] ы » - E - ~~ А А Фөш--к е --ж- ~ p ії s 
А oie ee м nU ee --.. мэт ~ E v е - е —- = =» = ~ Є ~ = =» - ae ой се E P= Ён - m ae eee --- Эмее-а pup T 
-— T - Ч мот қ. == Pu E | ~ --- - n E ES ~ » PES m М = аң - .... я e - - Ы m M aa ... РЭГ" 
E mid РЧР Тэ... ы тий - - ХК = - ind и ы = эн ын ^ ын вт “. - йн . ыг оне ~ - Гай... — м є, pn iul adii Р... 
ме з я -о. % -— —-——— ^ - - чуь ~ 2 - «= ae P rS Ф ^ - Ы » РС zl "E -. B > ~ . 2.” ғыз о ae mites ae --ө- т 241% 
“~ v ə ~“ - я ны аж = TI .. * == NUIT бегі - - - % E - ^ - , > - т т 2 а мыйы . 4с» - ЕГЕ - LI Жэ. гэ.” і - ~ PR NEP ntn ШШ, 
am er. mm mm odo лэ 48 ~ о m ~ == .. bud о ~ ээ П ^. 5 Ы е < я. Ф + = е ы - - й ер - - бу. ee і оцін ЙО r = 
рр ЧИ „Пен ILE ъ - үү ы “= ба > = «я so A a ә - . М “= Lais Б ы = 49 = м dile - СР = » .- -ж-..-.ә- “е РЕД Алы сені МЕ 
„ &%«д о ru ~ = ПР ы ps - = w= ~ ^ - “ Р ~ e^ - рг 7 LY - С . = СЕ Г , # + re dias ... че РУ PET T є npa am € arma m mom - - --..- 
"A P мж” мэ ~ ~ « ЫШ -- ЕЈ ~ > = ~ ~ Эн ы. М т т => ОЉА ~ E т е Ы ” +g se УУ“ КО. a Pes Ле Фу» даран = ы - = =ъ 
эн. ны n d p з «ду ae ПУЧ = = ~ • ~ ы E ы ~ Б “- нэ > » Га =” сэт в. . - - > + ~ ~ о ” PII I У P Р «тож лс - БУР”? 
ад е ы ЖЫЙ ашлы У Уа ES УУ мм a - 127 E ` ь L3 ~ ` IT ES ^ . 2%, Б E - = ы т P Ы = - мл PS ~ эж == ==ф= E - о о њи“ = P 
арар а мж e Ко ы мж... Е = э ть ы € —-- ын - ын - - ^ - р e o - өза ore а , - ы „ы " wm mam ove st e. P ПР = vos * ~ ” -с. 
«Pn зам gm Ыы У ~ ъ #%ъ У PE XII des ra мл ЦЭГ" ~ ~~ = " E - т Га ” “4%” т. ~. 28 = ол "эг... - - өл ә ә р о о іа ойн ре руча pM УДА и ~ 
—— à Р мн, ҸҸ РУДАР Уч, Т с ч. Aa ë Sri -. ~ си N - | ~ ж ” ть E Led - La — Ы = — ы - =» - * r лты TE У 1 Леле » on hee - -— а» 
uen" LL be 2 mp a Pipam tm - .... 55 а.ж. эь ч РЧ : » - а d ы М one ° ` - or « ӨР” ЈЕ E АЫ, ЖЫЛЫ киы АТАНЫ dla ји рт» — = =. 
aS "b en db ПР Ка а ~ ~ О et ~ + ~ - - e ^ + E и Le е 2 шаг M ы ~ x ма, Аррана У“ По заћи a omm «ай ре 
үрчи ее БР ь ть як А а „= =” - wA MM м Ж .. bd - 4^ - « = ~ 4 E ж - Га тт у» ы л С уп ~ Li Ке Рем РУ ХРТ " ка = atia E БЭ РРА Et йн n 
Mn эхн Y мэ рге с Ба лын tæ - ^ A У XI wAn ај E м La * ” ма ^ saillie: е тга adit E. ove Д --.ж-ғ % РАНИ Ке йке айны p" Mund 
«попри ар ND LL аны nase ря ОР ам р ~ ~ <> ШЫ -— ^ ~ lies ~ . ~ = ~~~ G ја 3 У Ld = - ні a ЭХ "ДИ т РИК Но ці — К ed 
= „= кирчү чү "Жы ~“ wá РО Чоь адд ET ы ч, ~ XXI nd ~ - о»? Е ач L ~ EM Р . Ф- “ЭРЭ Р Ра М т La * Pd У P -- РР 22 aro tt oe que. нів rs nm n "Il Amo m^ "У ET 
—— ~ -- n n P мэ Чу. „е B MTS PIE У ем 2 - - ` - = ~ „е E ы -- « И и «a = РЕ 0 77 кар ы эн muon те о ғ” 
кенет ье emo m m жағ”. Ај “ EET ас о ыы маљ А - e Б e ~ = е - р 27.4 вач x ” - ~ не ui ,/>.--ж 11 ~ > = ЏИ У и a ДУ жола pu 
aly ar I ЧИНОИ Хоја срца n mam omm о. ч - — анд ын Е an Bere — РИ Soe) * ын A -- oe -. нэ. ода MM PIE Да даа Pru. E РЭГ, дня potes 
mpi o = ы ee 0 ee Gyn ee TR БЭ “ТОСТ охь "Жо ^» o^ ж ы ~ » г wa ~ . -. ~ nid Ed ” po P => „>=.“ с. г - РТ шаг. й - > ~ Жабр иу ану ШЫ ИНИН ЖОР Жоанн Р Pn НИ 
карау Dd pues Рама а ыды r у > М E PT «ч E - -. ” ` .^ . ы wea. did 4 - цайн ал ^^ йана Ки - Ф ШР УК ШЕ" БЫ И a ДР Ро 
— epee try Pe ST tales ПРИЈА али " ч мэх жээ... hy Ay n ^ Та - ы A sm Ы А а даће че е саты же саг 2-7 ТК. ре ЕШ Бы, Пи РА Касљавћим ctt ira pm 
т а LE БУХ Жы а ~" ~ ~ мг” ХУ eed „> “ У РОК Жы ~ ь ~ У ~ = 12 ” P .. PI á ad nd aic. фета РО РО udi ү Г -.. “ сь рію» тм УУ Р РУ_ же- аб жү мем "m. - A - 
муч РА ЈУ“ ЛК а ПР 3 ны тт А о аць ao нэ Та • Деме - .. ы A ~ - ~ ы Le d ” Деба „б см бы 4. ке ё РР 24% РЕГ m pope nih dd РУДЕ ed mds am airs ei apt 
sd! РАНИ ЧОЗРУ im ртэл E “ГТ. ыды “.. ъ = Lu ПР iie = . > »-а t ~ "^ » ” ^ 2 ај а...” — io - e. imd кү 2 піні оба РУКА m peu КЕ ККК eee P E Р .. ки ee ae ee 
Иши ае рч ете ЖЕ, РА Т КГ Жырды ТАС ЭЪ, - о: ње Тэ - ~ t E - “ Td = л -- Е - =» он х орі ліні біз -— ДР ДР ae hoe у “АЕ «mtm dat t me ан 
РОНА За ара р ЦД УК Борну coe re Pas а == а ~ EI = qma % we 27 д а. е = a d ri eds MICI eed лыы ыы рер» уус" ЛЕК са фен РТ т Бы жі ” 
EP d ИШЕ" А банан РЕЧ ° me Qo ah om * bol: ` > ~ ~ - ^ нь! 5 А 1 ы? heh LS a: nm = езі РУТ И НИЕ У = ан с ада 
ИР а» МЕР, рони Фоо" EN мМ P abt: ~" T гьь . z ч - ұз .. 27 ЕХ: » Те + от e iod. Кылыр n = по d Pda ни up У 
МР РЦ Уве жал ә rra cm md m rat В IO С e^ tet e <= ж «Чт ын раға СА. Кө a PrP’. т. Фо чт Ты „р PE diei n dime н ЕН са Жа р ЕТ с 
Ар Иља ну и OMNE С 01 еее > - м E У T CT » d "um s ” = = = d eni = 5 ы ..... ПРИ јачи as epp cT ноев РОДА КА ан 
p Ж. рч, semini e TE he 8 22, Фа, ~ є nes Ч => З mae = Ц ЗА Та ~ • ~ - Е = С Ати. smp: rs -т adi Та d ПР РРА ак РТА у НД и РЦ Уа и => 
РРО а И ИМ РРА У ы Т, “ж. ~ mm. -> ~ о" ы ње» - Е ы . | ‚ м - а я 8 ee асв зай qe tomo m m Г на 23 БРЛРРГТ” а imita үну о т у ү =, 
Ранее е РЕ иж аратар А. У ак ий АЖ... о ба a ~ Ы Ы hd ~ & % > БЭ Га Md чх ~ E р ~ , . e 4 zl АР аон E А Да үү Т К da miii yt m a У бави - = 
- 14. D чо ЧК ke Пра ПР элээ “ГК. - - v v T T - - K а ~ ^ ы ~ ^. ғ = е ша ПРИ У ү Чү Т" С ЖЕЕ ш РУТ ок РРР Ун ИРИНЕ ntl РР “= И и 
w P "СМ ж. ж шқ, 4 а” 4л дарива сати Мо т У Tes дын <` КУ ы LES ы n а ч ` E La - oe ы Ы ЈА ... РЈ Доћи des: mtu» a -- B s eb mili м - Рр АРД К с „сй — 
ee е При ВУ ЕЛҮ К "гам i. P tm tm iy би Фе ~ ы “ХХ "be er On mS m о) ~ ^ гн D T ~ г ^ E - = а ыы ы аи v. ma um Доба даб (tse remet M T Pd - ча? — ці ON tee РО Да ћи 
а са. бес МО нач ан ФЭН ПИ У Зи ПИ ЧТ Тј) f “ТЭ. E ~ ~ "I «Р - 2 E ж "и > ы ыг i " Ра a. Р - ” пе е - я тэээт» “= И]. мээн udi РРА E poda PET ee ава ow ge a a Pd pru mn 
каға AA nm мада МОЕ "Р ы eee s) МРАЗ“ Же Ыл” ~ d» Ф № = e - P. Са n >a 2 =. е ғ ы Ра Ры S т „ч. кирүү а E | и миру у Же Уе = т бы 
en Ж Т дей ПР ЧЕ + рер У лана tad dd АД etm <= ЕЧ рр А А E LN ъ= е6 - ^ Ld - ж PLE К О а ее ына ЕО Та a У Ро ш ПРИЧИ и учар ЭР a 
ПР элээ P ји а АЊА РЦ УН “= Ww TEL жи + + чт үү „йш ~" 5. О ЧЕ е LE ^^ ог 2 E е --т ”, y г 2272”: ғ й РУКА. МУ Т К ЖЫ жж ~ " РРА КИН қа У КА Коса оће РУАНУ ин ИМУ “К ee E a 
-—-— — X e РР Е "ы. „ы мэ UR А E XII eee A нм ць wt ~ РАТУ" E * ode 78 кт ү КЕС з Lad om ~ ae ` * „А м ^u • “- ^ и" эн “. per ТРУД Табан ы 4 "e P сонні лани d тир ой ү ЖА мар И Со ЫЫ ш 
ЯГ ны АЕО м Мр УУ о РРА ЧУ Ч од — 2 тк A ә со va ы ~ — e . - С. G E or “es * pde robs КҮ УСЫ. „ААЛ ПРИ РРА У У није „РТУ "У ТИ „ Сы тнл р г б = Се я 
— MM ern aa Яо лэн a a a L PY RENI Ж ~ - A те = - 2 we LI №5 а- тайлы pa P Fas a ae at rem t indi nal PI У РК aedi РАСУ РК, чүү Ж. Б ПР іле anata ЏИ њи - 
2 мэл -- аа | ЕУ о Жы ы m --. - - 82 - 21 = т» Са av 2 -.. ғ БСН, хайлан мн. о а ТИН ДІ.” aaa ЕЦ ВАРВА ме = + ун ДАН план РОУ Ги 
шш камен ҮЕ ви ме” 5 ^ om e Б ^c UN А. ь 21 «т Ч - та У AF en ач СК СУ а үү теры УТ Жы >= T a Та =. № У и а и Ре а | a 
~" = У Уни sn ie Aq У in iaa cite Lum ~ - RI = Уље 7 "PT mad ун А ай жа ww 9.09 we ү ү уг 4 253 йн 123 РЕ "Ти лын ни at рано m БЕ аы рр че 
: лер чүчү чү чуч "КА Р А" Дыл Чи ње -. <%<-4 ын мэл a — y .-” du ~ ль © * г У, ы ө" ” s б ДГ. ГЖ. = урин У ДОЗА А бањи БЭРОКРФРЛНРГ ТАРАТ ТА, Дал еліні "ар. пор Жс С Жм кууш ү т дген жт но ба 9 а 
ч © aan У te a шэн Ч, Se wee, SH EE Ele t + 2d i - Мы дейі „© „+ 5 “. == s “~ uit. а - жтағее РО АИ КЊ менее” Буы мі ^ И жі ча ч= а екенін үтү ай albae M5 е чу 7 
УЧЕ терчү чүчү С ОМА disi чач ч, "ЕТШЕ oy, a yy Р e Latine Pets Қ 2%ғ% У ~ ы 2d р В E - Р "T 7. Pd - * АГА Ки тео поз ба я ж вро м Цан гэгч In ames cium tun m ма аР Р и рар рр lg sm ade 
à => 5 ЕСТЕН ауен „та Ры ч е У << 9 зи - ~ ~ ” m ы М о > t о 2 Си A "Ба КУ ЕУ. мэт ни - ПР івано ПРУ У бана 
ЕСЭ ЦРТА РА У У ОЉА ЊЕ ВИА И О и тэ? TOP V pl asl ъ= " % ч 2 E a A ” ПР Ы « d v М p dms ршн КТУ ЖЕ Жыйыр ee дели РУМ а о и пошту мед =" pour asse dom Жыр 
dd — муу” ы М ы РОТ 4 P T р ч“ - ~ mec AIL . £ y о “ А ні мій майні шэн МУРЕ У ИО ЛК КЕ алаг РТРС А таи | РКО, уул лт ү Лс ОРН ДА Рай» айын 
5 5 —_ CO ba iocos enc аы. 5 d dis x © diy re es em eM ML a Уы р а оа 5 ғы) B Р a мж E go tmn ET Ф о айы үүнү ы у Кү idet m" è orno qm B am muro geram A m Popma pena н алта „Жы 
4 ~ ДР а ан и р парі cta ug ЗРАК ТА ын У «n = " ” мн, E E "ricas Е г * ” Е лаана ===" МР Аја ө о “' = On quoter tata рее ава 
ль цан > И од ль лы рле" об d — qe cm й м ње Р РА pv dio puppy p debui 4 БОРЦИ adsequi aan диодын И ово 
Р ~ gq Cc << №. > — Па нь id ==. КАЗО чо ы в = € Ч К ж. О » ls ы ам м " ... ГЭРЭГ: - > м СК) РТТ ХХ ГОХ... ээ И К meli quiu жем күтү ЕТ К Шу о Жем = "ЭУ Уу Еш СР? 
= ый ЖЕП ы ыы дады " M—À үк" — ЧЧ У В Pv ^ b EI а. 2. eaa >” - а РЫ - РЕ аланда тай Е - үү сы ЭРЭ” дэн -—- үрчлээ ue. abt s У did РУАН ien cell ci 
~ о = < <Г ЗР с ы ТА o ae A, ber То сб йы даро оъ +» Зла, ule а У ты ГР," Г, б М түй, Гана -.Х...- PT Дооел тачи РАНА Уб и deaf wh feug o ab o re gue ie eel ae eT рат керн 
2 шоо ran ааа 4h 9 mdp dw пој eee! Мас а, = о ~ s — ын .. 5 „~ а ны Те P а ы Р БЭР" dinis Ы є У a9 . таса мина m ut 2-27... 22 --” os туеш" Л К К С СЕ д» ч© д” (Б%фе* эш ан пиг ал tie РЧ И б de „ 
~ c = СС А т X te at ee ы LII л УУ. о a 2“ due ^ - ine , „ - P ai da n i “а 2» ићи РЈ У Два АИ ПРИ ТЭХ хилийг pep rr кеті „РА њи а » 
TEM о а а ү ы ‚туты "ТС" da ПР ыла цэн очих е = м t% + т AE. P Га rey .. эг» sesta БР КУ Ыы КАМ NOE P an ud n en иш un ead ne дина и p 
һа Е LO СО рыр А Хоф p bit м ds ined Y + Е т ` a =, 23 2 2 М Дай» а гРЧИНд toe РРА УУ Гоча РРА Ревана РИС БЭРТ” 
т. 42 — pq љиљани зэмлэл” ји ~ б у 2 22 ^ «= 2 ы ть tnm ne = v ы we ТН жағ ғаған марбырыар И барчо ыты dtl d peer Жм ПР обр оды 
“ a (О і пар tt t - "peg, ЧЕ Ae, Ca rt К С“ LE. >“ T © Ф% әм P ~ + ~ - о 4 Ы ” У, а 24 E УЕ абиб РРА ү нь ни МДР уе > м Ран и дим — o! 
— РОЗА Ри a Y И У њи "ы з ~ • * ~ * М .... = а ” М я adii es ы аи C EN aa абчаі ог ИУ ийн алийн ыы Ipsis mad тоф зу и mco цінні 
= D ENT У ен Ра к ов т е дылы да А. мг .. мы е, Е ~ v 7 й Е + ^ Lo . 9.9% o РР ОК кий Ри perna aiiud pend ap слао ваља сања 
= 125 x te fa дн еліп nme PI реве чета РРО ОО eh oe а 4, чу = - - ж у P Tl г. М amr ~ “т Ён "s ПРИЧУ У со очај рф еб ба им ани үшүү" d m 
ages = © МУ mi pu acp А, 5252, ћи MUERE eniti el Те Р ~ + 2 ~ > v > “Јо ” я Fat dT РО Йа к ЖЕ ied OP ү т к ш РА Го. gti ieu Dt oft uf sa MI рељефи ај = 
Р” с о "NT ы чт № > er e дед» odii adt Mr а а ита пиће wA ам ~~~ TA ^ 7” д 4 а УР” ын , rere a ees " LI ner dnd БРР””Э?”” Хэнийг РРО ЭЭР” Та лаа СТО ene puma а Ри Пашин 
Ed qa © = 759 Miete isle РО i ~ Фоо е А7 A. - тт ы. 22 м P ^ = Luv d БЫ „ ti da Га е. т d Га ae А Доња ићи Par Цай ad үүтүн к жы ыйа егете ке РРА" << La nelli ed жие њи мм элий EL 
ae > (түмт ed pela Se e ains itd РН y P a RE ..... уіл E е о“ A У ж " * irs E s ^ л А ..... E „ avr Fe a ~ n а 0ч Ж. әй імен я» ПРУ У «и ће ГҮҮ 14 13-44 РЭГ 2 2 аи он «ілінеді putetis и мой 1 te c 
а © ““ 9 1 em к Шойын ПРИ Дери — ob jn чить om oo B 30 мж. №, - ” ғ“ ~ 2% , — У 23 Poss rM a "е быра ыйы "Р к. чб — ‚ЁЁ Кы ym suntin Поинту mo atm o gti 
mi КЕ И ви eee К ы ое эь E қ à. - "o ем й — үл ” ГЭЭ А, - go se afe ig rm din In LP AL РРА Жы ial РАДАРИ Y adiacet шщ а чо оч Р 
— — АТР да ды ЕА РРР tinte in Г........ ---... ды ын Це | de LIS 9 ЈУ ^o» E Фе e Purta Wada P PF E .” "ig S Pv -ñ РУТ ge tora Р АЙ — "А Сања о Ды P reas Or derent ce ы d f c ПИ А qtd 
СЕ со > РКО Та Де о eh een ad х. ad е чы = эь. РТ АЉ Ае. = s ы” p E м ee "Tit * .. - pnr en ж СУ RT К е о мо діа уч «му ТЕХ Ми УЗ жау Дар. КИИИН ЕЛЕ Т Сау о ы лен нені a ond 
X L [ і „Фын РМ ~ айай үтүрү тт. "Ыы с >. СК а EE dia) ПРУ У ај ч Е * E Li ~ p e йы P add , Тыа даљ БҰР Т ЖАРА йр o Ка ПРИ інь ліві лай ее ш ремет же нений БАҒ Та берч pnt 
(o (О де љђинији 4-....-...2- РЦ У Но о њи У њи LET „~ а (c ~ PT P РЕР d LEE Mad DE ЖИР ТЕД э р ЕУ Р с, qm er ди риРтт””э””7” m МТ ем р тээ" д сей Ар ЧОТО» РІ» екін i иф ami o а татым уын 
"T — ] Mc e а аю У Цан зарға "2,2%, Жы рд ч Ц Мн ~ — ы "= # P rau - т ..... ГЭЛЭГ эс - РГ се IS n пи I Fo а нні n d car apt p om at qu aa at. сн ИРИГЕН ЧЕНЕ = с шуу == = Бы Е 
- Се ГУ Tcu Ф. ето баео ари О ЗА ЗА oe aS T hes ~ a ћу мя ~ ч e * Р ы С". в Pes w Pre ~ [= Куе а СТ. ИТЕРУ ТАҒУ т шісі б мерт Т а МР ы. т” гай Ч. ай» очи Pd инь анжийн 
- су Фф = C ДИРИ “= А pP фор a didis e ин - “еу” 2“ Ур“ У Р РУ, — И i A - Е ҮГ? ТУЛЖ оч pep qn с ЫК б eiut il пыльца 
ай, CN —— CN =) ben pense ы МР ni d oodd — % “У 1. Pr EA A A Б" а f > A ғ "p ” ~ a 2 ЭРЭ” Zara Р P "m a^ po t oe РК РИТУ У а рн СК ћи т ғат побарани е ігі Se ae Берген == 
+ ч Е M Pci iie dn МИ оће мэтэд ИРЛЭЭ Те ы Ae ТЕ. жд. , ..% ж- = . ы мы 4ъ э ” ш шаа аы т еле ы „= dil „быа y cou е o d n gnum ch ue ames a e ul at ae m e У меко РР фино ү ын нв 222 
8, Жы! с... C К. a gl, 4Ҙе. pu Led и ћи к” i ay RA 1.22 оф а II мы. ~" ГАР. М че AA Pa A LET, m 7 тд А, ж” л Са Й e UL abd еі "edi ditiis ILL Гаа db а една о о арен arid p mmm mdi 
ad > ЗА. ты чик лин РА ји Ён аль boi атын С: “= СТ ween УУ ~ “. "TP а У 2 ш Р А» айлағы "E здає Pe Рова РА pod and au Раја dsl cn ийлнн РЦ Ри e^ ama o parque и црни У н 
~ (27 а тш e a de айы да pet no Bero eiii X о S 5 Маји. ` " re v wt Ww “= өл "Р Фот T #/ ж ж МНЕ, йа! ређе“ "ве" куу үү Кома »”” 7 ыбы Pa "РУ ак дели <= id дар» літ ГО iai if porq амилан 
М o pretii aL И Wes "m За КҮ бы t er adf o do ut P, P цена TE. з A LI T— я " Р Poesias PT ET dd А, У" БТС анин Par di iind 2 и р иг дэнс A о аеро md 
Q Q. a Ң/ wt wA У Зои + аы ий маг” ЧАГ я L ~ тэ irte f УК ” .. х A А 2L) pa gm uo ame wa P A "aera dites adi о и М Mb sl ttd дерү унн чеш” „ы ринро rim Rm ПРА "Сы 
2 х ji кізбірдер қулана далы кі дін ды Мар ер Сиатл tort wae M WE АЈА а ~ > E 2? г» „“ Th 2 о ч .... 4 rs EL тэ нд Pe ТЫЙ ЛЛ УРТ К inte Ен бини Мом фри рано Ири pas ни 
4 = mmn ГЕР ПРЕ о ан КУ ал о мэ Гон. ман... "Ё ыы * . ” є n al d УМ О B EL е. ы У РОГ PI "Ти РЕЗУС) наз "Our аьа P ii i d M Lm mimi cd 27 07 ПР ба ИЉА и 
Ј LLJ M б пра и і dicam рт 2. ын прућа = атъ т А ы П чај LI ^ > P У“ = E э e 2 А я a a a а .... А а буй д РА biais cd gti tati dveyh stud quf «19 ue a n map an P РСР Тоа ладна аланнан метан ана У Г-__-= 
F. и жа желі ы qe нын салла арда ды ee [чие рүү е кы. “.-. “” ЖЭТ. ын ж ыы a №24 “~ w A тог тэг” » <, РРА РЈУ адна "AL: ne? teh refer = See Aa дин итрај и фијина Meu Ст AP = 
- аы меи амыр по А тн а ініні ОР, ” 5122 уз Гв ” ЫС 2 ман o - й маг”, T "ENTUM нні - У АЕ L БА d adi ем а re adi слање раро абе ut abu и ДЦ и ГР цан и Perd 
d s УК мече А рм Таса. Бі а мала ры дылы ы „о ~ з № - зх ы р qe SS Ai “..-” МА Pos ГЭХ e ь ан PT m мот ДУ нанні най РИ” ы санан алты оријентир ВИ ру чани pipa 
Це беке ака tiia А ој РУ Puro e дт 1 ње: = ы BN 1/74 ~ га. . P HEU d = "T М т үөнү е че" ЖО К Аад teal рур УЕ айдала t Рини РА РТА аи а одао. дын у э дий 
КИРИР УУЧ альн UP Тен а оное А ји “~ у R “У 5 тээ... n У а + eg rw s IP dai ог РА УК ідіом siii 5 күүчү d і S РРА У ајд р quai d iod Arama 
rM ће њи а ose de rag Ort rar er t hast eS ae ЧЗ. sc. dv den D АЈ ад ~ ~ m 4 ~ ш- - Ра a rn " Р мтог-эмл У јан ччыхЗае 47” МИР Т seh ol deh л айий ap ey v gno go ЖАР а PP Го РЕДНИ У УУ об кичү ш ч Да 
малим ара Зи ас a A dk y й «ы ЧАЙ АА. Өр ир НИ ЧУЧ та ЧЫ „ЖЫ РО РИД У Эдийн И ~ am geis Pas р DRE шә ato y ү, У Ка ем Да Ta T ai ianiai ж Para pni eter el inii пр рана о учора у нам 
У NAI РАНО РАН њим тн d даамы “аль —L қ ҢЫ ааа, ыы» " ц m 0 м. Е. ти pa A ad > — „ к К ТИ о по" У и а али РТС Рм ВР и painter ta oie rN aia — edi ` 
р ьо НИ ји etit оче С Ари nad К "deb "nae У ТАЙ о ~ << бы - Gas » + г Та -- л эй. ur es som З УМРЕ РИХ о о о нари амур ато» сан РО раро долбар раде фон фо иу ані іно 
n іно rn іні реа at ићи pei ам ДО Ања ч ЖЕ т у ы ы, п ымы; wA c o т w^ d 9 қ ы. Р у D = an hr ПР а РР ад же ім м, алланы СЕЧИ РОР lad ale eth dle ыы рдү" Ус ~ а = прно У сат s 
иа МИРЕ ЕЕ “ы ТЕГЕН "ЧН е А АЫ па аде. ЕТ Ма им Pa ы Bes СК е TU M — РО є УУ "У n ои РА " fA pardo: то є Р) У РТ mar РОГ п ae heel РМ осам "~ ПРУ Раца Г дох и лид он onl AW SS 
ће nnn eh ені а И Таја ОО анынын eel LEM. p wd ode e е і Е Ч "— "Lr ч "У wo" У мн «РУДУ СГ У аи Берген о Р і вніс te мање od А МИР на ПРУ ни ср лири ы рр 
iiti acie etie ны Қа. чі км CE a и им җы ынны гі ар Арьен m 27 Pee е агч « $ РУ Им К, Речи Фе o К P эж Р IT P niis күни ЧҮ УЖ ы керилүү ей фрун ө Й petiere prd жәе ppp citius 
tr обнимет РЕ М „м дуван б, га сл se на d dua Tid Eu а МУ eles ets “а ЕУЕН ы ы, дулый a р лбу < v при и di И Фм € ша а Е. - 25 smt an n Tn - upra fe fu f ml PLN ы тт на go У ји >= жама Ра ада dia dtm tito И а al s, o. atm peu dimmi 
nidis price an debited: ди ар ды i I а о ^g ч- ha uias onu, Ap, ће а у Реч Tul! ж-- “ «м Tan 2 eum 7% ГА?, е 4 Li m “ЭЭГ Л... À Т і на РУТ омање РУМИ УУ и МР ои ОИ аа ратове пица денара чина ма 
ng амыр ЧИ а дионица vw un wie е сите». Ми а ока њи У » 4 а ум з "а а ыг . I "У Я го, PE У у Яг” m = йыр ~ PT У о ана РЧР о мич =r at Pre јр РАИ eaa ara Фараб D a Ea наредна 
о aa iA Lg UR n cA ORG a ATAK реа ен рсете А әлеті» но о МУ. Мен дима р е и мг, E " тээ” г ё cm oe p " ps Hm X wey * тэн ай алаг Жол өкінеді A Рај Деда роу ај ти ари дни. ПТ ДЕ лан 
Мо оф poe ьо T аайын э а. a aaae e 27 ПАРУ ЧУК Зе tA тэ 2 -- " , МОРУ ја. 4% Г. зог “ md mE ГА - s СУ. "аб p уме оо аи иран n 
aiaa airne aae a SO iraia aR р i r dr 45 aa ын Ж ғ м-гы үү" "л и ЕЗ "fe Ñ к Von ч тэл Булэг" рт" ...... "Ts E NL. агро чт бо ям ро ВСУ мй ПР ДАТИ Ела а Ка о, i ш ыб ен НЫ etr Qui api ром 
Пир О аладын |ИМЕ іль Jap pta, s - = НЭГ” ДИ a adaa i > ы му Ф зок ФФУ |. T “и ... ~ РРА АРА Жы „Аад PrE di Ки Фаїв Ке т айбын корүү чуут" Жы ое diti un ычча der 
Е ан дан enu E И оный pa at irrar лови р РАКМЕТ" Ж бабы ы м А, um о Лаа =. з РТ a 4 # od сега Ф! . үл 1. мрена вео дна Аб авд Р wot шыны тыра К сш финале 
О sant pata etl re —— ibl њи ЛЧ К ditte deine ты” жыды ТӨГУ, matin ie b 2 v WE "= “/ "m to.» “ ^ тү, Р "LT f» suff vt ES tao t Тэ МР ТРИ deleted aie РЕБ ады 7 Гани ОТ чараен: ыб И cirea 
а оваа РАИ ОУ и дЫ па tea adito die ride У и "p ам TT. we г E . - os 17 "anas. a е ТА nah рер тем тте де ад ње чє а-ай а на И ла мінні ні во р си радаре јава 
ПТ с a p инаћиим а а А Ио пр р ss hit in кодів У И ~ ре арыла 2 2 E УУ ~ ~ о „и P i: - ... Т8 PLI PL e мэн Хн. алийн Алдар уыл њи ат 
и р нар д билээдээ мэйл Цин Эн ји Зи hg te een emery NY Se РА РИЧИ У ЗА ОКА ДР РА УУ ом РАР У У ~ бољи 7 р Р , - D "LL nd ur tete ПРИРО ПР К А жебе nu и сања азимуту E ооб узан о ооо а 
правили Маја У пад aiaa Про Фу бо А мета об игл 7 ы - d T аа рт, порт По чыў РУЗ. у ү aed к” ка Брен РРА Сен РУТ Ако е айбар САУА” "аи bud о м1 ра рр араа фро tait тории 
а је А-ай б дыре Mites мре ЭЖ, мг 22 E Р оф — «Pro at — M " - A Дуо 1 w w ыз 7” й . мә Р, "rl Же ad дығы, Үү Ald POF A OGL PM ET | РМ А dint айн н порно SID. 
parr i гэг өл ий а а о роеви е-е ође МУ no, e adi e La a m AL a “де - hee“ fume о 2 ooh Gen © ee - о. ли быа дй а нт рада, дай ой Ра хилийн К НЕЧ ы а айный mnl ptite direi dii gum pet AITO 
erp date tt rdi pcr e t ending ier A „® + rri erred ЕЕ У" Эн. ерер P at ou Pa а UN ы 3 A LP „/ T T "АГ. ү ыы д ~ "БН ad mr uni ОНЕР И Раље СА залны ПРУ ШЫ ДРУ РУТ фа P = Же жый ag arate erm Н ныд 
peste pub grt анна | d Mo Є mache E db bi dod ГРЕХ __У~ ep Pip it aat c.l t, vds e T андын дын 7, NS bra « E = “ Гүлү". ar РЕР СК Деш Ч” и мг а ЈЕ а ПТ o m piod е ы Жии реферише 
хань ца дэл дуд-Д Ан Али Эра a tah eect ere tan eh ey doa Ss oth аА үүчү aa Soe hal anette рэн E МА ТАЛА а Л У У чел У V "d — P мы Ган. Ра Dd У рова ЈУНИ РЕЈ ОЈ оваа IE д4 О omit о ei bulo piti ч. 
Tarara apn ыы, ћир дна ај У ная Nee Ld муње њи да "У Ур "У “ РР he a "L^ 7% „У " ре - bas i n ред мы Т но А ды ilii Abiit p ft tiep ape ренина таван 
До с и m Fan dr od ur ей рм nodi eta РА и rel 2c ah. E Lii di ^q Uy y nw ~ ЗАДА ЧУ А ia лы "и e "X ? 7 wat ELI - P ” ” ИДР didi idi мила МР Трава т 
Я цьо пари 4 pir RTE ew ye WE mi rcu ear. pP d ЊЕ У S dinis 2,4-7. У терге 4 уље" pu MP РН т Tum = ә, Фе ы А Рі РРР ьм РРС А асал РОНА мд К де о айн i ti оомо cp e m 
——— ль lioe penes edt iari АА а лу до дуйрйн доно bind E abd "= "fs So rrr uh Rn mad РК У "A TE T ET VIE АДА ~ С И ЭРТ че? Ps меті” Р " до а ж өз ҒТА А жа обедов бі РРА ама Уа реноме tiran рефрена пива ај лана і норе нове Б 
Meier in dt К т ои Мо ли mou бы Поља лэ Y PA a нь О ње“ чага пе мас ` T Жа ы бал " PN У ЛУ ү d ы ГЭ? ДХ тс Ро Дана Pd Ра Рм ФРИ eoo € I и pisi dini CO qr factae toad 
јами no ing ri Eata a iei AR AR er м қар ч Р аж амны e нэ Кы Йыр» фы v s Ма А е. Р нг зан ra си РАИ Ал | qp aan Жерш m PT d iesu dad ПОТТА idit Жылы офф qtd ir nn y vn tque i titii 
о төр ыыы ные Г КЛ нача ЗОРУ" У“ “<= Ру у 2-а ? ба d > е К r PE E Esai ted Lati b a дачи ико“ ар ыйыы вала рана 
Чек пир итал» қыр се сан УНА А дођи аф P DOR esa Papam m pe Поа eed ми азия w^ Ж. + че У m , XI У. ч М РУТ УЧЕ ла авы РАД id ini pai pretio етар пр она pihini pagara qu E ли Б алд 
nrc rui tiii а дачну a = Ма РА ы Po ам МИ ані e nen eC diss i act T€ 7 ТЫК» ымға я е. 2 -» [ч ТҮ: m 2 r ПИ је TE ut M d Пи МОРУ У наћи ея 
ее ite крш" ам Me еф ан aire Қ ТЫ ыы. Жас дын t Pà Бо и о у a УУ Олон , 4 ес ” ме w т оба РУДО eil ць й ПЕРРО ад дні дов Й Га міц ікрі очи «ор e gabe iios у Даљим 
Я фа Ал Ры ago Qiiod dir У s w^ i г a об Р У“ й Т6 ee И, | РК ТУР реа Маји ај аи PP ЫР өтемеу “өрді oom n di d јр овисна ране 
tired Pa ділі А ји iioi in FM Жылады, 2,0-4. ic ларви жиг адын! m tam == Ч € Bl a aT, “‹ ' « СУ УУ рам офи ар ее дене ТА, ipee icis бурд о ық БӘГЕ чыч өч йз ра {Еа ИИИНИН НЕН НЫЛА "дый 
Айды ча УЧ ш ы (6 "ws У њи А ОРЕ кімі а көй төзімі e МБУ ач + ~ У p ғы е , р “У 1 "гоо? р Ф DER. ні "PU РК аи Пали ET D ETE Tree EE љета тег 
рі з чинні і he tel кю ані сф РАВ Уч л р д Xi За ~ amo ы о Де Р У * ла мэ ЭГ” ї е” с ax МА IET a "d 2 меу? ићи РСК цім а adel PP ON EE E вый peo Dandy 
nd фое оь АА ddr dudas lee eei ESO а MIB и dite qr nodos iiia ај над раје er n E р У“ т ығ» - a Та d <= ro sy ele ow мэ “ МР ә т согд тэр хил дэн Кини абаненту ода gr t банде срања абран m 
лінь З ђин ЈЕ Тода Ыы Ун" К АААЙ мађара ПР ње Ч У или "ju fef. fri p Rao AP I со 27 а А) 2 я РТТ ate УУ РМ РРА Дио УГ и pes ewe У Ft POT SOF n ub ud didus Pene ial ia] v eg oc Rel Qe ie tran келерін ж. 
нь нуун pipi he itt ів цею бен С нь ал | v wm ЭХ алд а а b d ть D TT EM ” 5 n E ч 0 гы + : Meg У ЈЕ аи Рани РРА a ATI рр aL Ar ial peat id a pageant rea Nera ee 
ee ap ar diat AERE prar eique rire A: ла о eR EO SA Pola P а у арт обед. Лү Уље Р “ П ЖҮ РИ РА | а? "ДР. ПР, УМ МР У i "uy d APF аи < LEE E ie ded ppt ip У а ш шок оре еру АР УЧЫ “М алыи АЁ ee ee aate disa 
ata ae рел oratii ae Gi EL ais bct eo edit pede pig decipi ыды Na a a. es nd hte iri ӨЛ ahah Ya Ж гс тебе Oa PA Ж Е У бич ti PP LE Lei mi МОМ Т ЫЫ Айн Е 
Боож Загаысче аг өэглллгдгчгэээц ч челүгачфчар таг эр ал ӨӨРТ а ди и ТА ur тања n талын АРВ ^g agi Кн са, ісі, 106, 4% пёр ~ . Цанын Сана d РАСА, "E ПРУ А-Я НЫЧ ү РЧ Доњи Е аы РЕ О ы П на мр урне а аена одолеете чудната 
и ји дыла бду a a a i ы ы Са rh ice dn Аа вн РЕ do wars ni Oum emn FS yea «7 ow t ТА о “ = eh S Pee att e Р А aie (GAP PPS m aan rom von Vea pit Mie eR ug id at iupra A Ap paesi atre rae YA t e d i Er 
да чьем predi ам РМ ићи фар р Че agi рр ААН лы МУ «& "onda M v^ P. тэ” АА We Л аи офа НАЯ езі” ћи Р бы: 
кс се рі нні hee и чт WD. име МЕ Еау "а y Moe Ou e EPIS TV "ыру де MIA LA dede 40% t ТУ. » t Ya Pu E «x Vl qeu pu од a peee тыч йс * cgit t А ча ГЭР ЛУ edi di i Bai p Er E STE EAA A ем ч 
forte фи оф И E TE А „мды еее рр о ИРЕГЕ ЧЫНЫ дэн wee О МИ С RET d ied M ч Sos њи à en Ioa f Ч?” т 4 фё БР „и aro cin аи РК аферу ње ПО ви fist depen rer али x vto a ш I ap c duas ag SAP рав 
Le aii Rei Ped fir i Ae ie le EPA B I e ЭЙ, ES obiit ii dod Y laudi P pde ancien ih end qd 5 м REV E аде "prd ode d a diluted РЦ b aiii ГАГ 17228240 ети Џон 
ор род tied ao My tore а aii inh Lb dd. ^a ferae te i.a Ye D bean ыға мч ы түс" Жы tas em yr о ды ф РА x » mnn ” АРМ E or А ИРА gh FP PIM gist = HAS, Мо д ды А асым ді ротонда Атини чы 
Ах ЕТ аи мм ари т "У У - "У б "= МУ АУТ о а ЕФ бун ЖЕ УЛЫ ЫН "н мі не едва о ротора финале 
ih indir ЫК рне E E d Lee а КЧ ЕТЕР À ДР, Р У ан И ИМ јин Го» а s MA T ds nibo з” Ад ар МА ЧУ Ри р y dues а м-н ДО фи арени пере 
о iA У м ме и = мб ати пає i) ddr dnd vt avo њи irt: ^ ira ote "ad мн” ЕА "e РР О ои ПРУ аад РЕР ао фам "я УЧИ reden осой файльо БЕЧ арр Ч С С рь чтоб 
аланы дур цн манаа лэ вена tk a tt л Та АЖ ч будл PPS f». Dew Pa e a че dere asa ST EE "Р ы + Im 3 Ba wA f £u o pre 4^ ed "T Ме ВУ Пе СД ем РА #4 fT eddie рик 
о rii e dide Nl heo i Nr AI eun dm AP em zv Ај УМ ge ыд rA а чи "E tl Лк. з wt o^ FAT EE I МА „е РИК" ү "үү. РР РО М С^ ЖА ЈРВ вані ива А ет 
парни ва Is Dada de rie NR ды пре i Lin A Ab ie didit pite Lyc ae, ЫН e o ам ье. = А і 5 Den P n nu gU V У, куер Н х бе Ров ари а байр р 
а ваља ipe edi res ene yr UD Ir m Pn Mode RE ede dur а туру 52 u^ T Ee iei ы ру ЖЕ: - ПР О Аароны P "dd E yu (um m CP РУТ ои edition вт уроч ма ад й ЧЕН ВР рвана нн а ышан) errare arm direi 
bag Co wy mre ^y fana og joie tini ee ad i dili d i ни rab Te д. eel Pond аф а е Та ^ 2. d д є iw airs Dr У ЈИ ки Pe OT a T te a ad dln с Р-Я Fg OV OPTS pa nm tte Ma dri rid i dei foot Mom уе. ж ететіне іні тәсі іу ыт Өрме 
қалы ала қ ста а Аы Тал АНДЫ ан МР Ола ды ди GET PIE =F n з мн E еи Й n» = Р zy vn M і АЙМ АУ КН dia io didi Mn QR Nw v lera ve ut y" ms Е 
арены Mu i ВЭ poesi deb ЛЬ РО ТУ a ~ фу Је МР Ре а че ау-ет Рај“ Ч БУГ ad чу f е ты pu ETAL RNP CF р БАҒЫ зона ere PTEE 21) и УУ К, ae Mrd aed rA eR riii ra get dA eri Pr Rae i 
ua a erui s P ag Myr A Ve PR P AIC a ыы МЕ МР и ALOT OV ч P cia Yad ho > ЧР оаа а-о Роа бка. UTE а ерер НАСА ПИРА ДА РТА Миа ден sea ILI поп јаку 
цифара мировања eit sth tt tale oN A dr rip P eA Ua M Riiie альо eS Ll rte at lin th ale Se менге ле” VG С, ^ Kv Ed БТ T Қоры оне МР Дени Р РИД йн | PE gilet rp a la e ie tA M 
ed AA e hinge d lca ARP M e Mr УРЛ а да M Rd о i das: we de и MUT г wher Ow ГЭЛ мыл” ағы ы өре» ме ee a^ "v рут Ус" Кеи, РАСУ у” к“ эч] ГАЙ" За. МОРИЦУ с Жыл „ш Ета 
а M алаг п Аа КоЛ, ыле ај a еей, Ты" дб, ыл, [Мр ыў " Фаг «ач т 1 т barr adag did yr ћу 59 ТЭРЭ одо УМРУ АЈ рат а памет ут рит adi ere t роја са 
diei de ae А А иным, и о чи a RII Le La Mag Md EE T E Улы тн ” Ае белум, буу fa à але дее ТИЕР А > рае м V бұғы ој aipa дивља ОҚЫР лас мр A SIT La a ggg ae A те 
ая И Р ар-р ы Аа үчү еи ы А TN a ҚС ыы » pr “ РН kd "Hx T ҮЙ А Ку КЕ ЕУ а РРА в m el м 
ОИ дары АС сч у А З о Балу £ Ји . + pòd ni УДА. , МЕР АРТ ТЕ дар e Pt очні POUR un ra rmt a ee deo LR ЕЛ АҺ, жады AU раван тата атое 
о einen bu) op sg e na Ава фи А а дао) а мм пре За дюн жолыма ДЕ алғ адм A. " ~ e " «44 r PP T Ре вмч als б Ls rn (paco А и ur Pa A LA di ыы АЛ Ат 
АН, Pee ai Redi edt МА и ма = ре АЯ Tee. ыле M. 4 ‹ РЭР 4 е uer var vio a raf iab Puls ta m RO NO ert маселе Ры ы PU e deaur prd eT И Дену ејд а m : РА prouti VT v 
TEJ а N T ла Medii" КЕ За А. МР Са эма у et ыг ы " P LT. m B > Ар е УУ. Раду у БЈР РРА ЈА А ih РЈ РИТУ oa rh n e rv ae = “ Н 
д аа Ф А eae jd ыыы Н "ЛЛ Ro dd ie и a . жэл, p AVE dbi di РА МР РЕД и. Војне И дана | а 
“т РИ ње А omm gy b gt ou y e^ a ч vro aw RUE ТЕТ d я PETI ae зч у = , s “.х... ҰТА Т һәмма ағыл б ғғ” че ph БЕРЕ Ал E 
ы 3 мн addebat Yn Y TA m LM К .: а ~ VY 8, я ЧЭ Ре мечу в. 2 (487 <” быт ды ПРУСИ IE ee дамын BBE CARAS WOO FEY 4-34 лы i Sta Poe OO ран 
rete и ићи ace Ww $ ж E SMS і "T T РРА ы e qutm dd intend Pr rn A edes РАДИ РА РРА дни МЕ дај р jade ai db tar коло de 
М f - ЧУДА У РУС Ан "d Р лан Ји "T odi GP di САҒАҒЫ М2 initio маты isl nd Орель й Бан 
Ы чя: ж Ce бүл УР Ут 1» ГЭ od а ped tte ah ұт ТИ хэл ane ҮЛЭЭЖ 
~ bero р „до Са Р но е ЫЕ 
on etii Ири на фин офф» 
льдины 
wo Padre te 


иу аи пада че Му ЧЕ а да 
њи АДА асс Ма а а Зеу а 
плита не мума прса ets Ле oar overt ај HY NN 
Асе И др дэн distor tenn titel tie Mitte) аа а ная y "un. fy ade А рее д 
ЖАЫ ЛКЫ у АЯ) Та а ји ла до ва фи и + V Р РА 
) а зас < ЧАК на ма Да нг, 21 =. 
А ы - % "sx, $ %ь FI E acad disci" - Р - 
1 у „= PEN up qu^ vp us gr m eror el жилд» 
7 “ИМЕК voy m er AV ue p "v breue n T d hood РА Рио ла 
РРА НОРА ДАО А хе POT Ты айсы te ae М are ПА Мар Еридан 
apt т e bendi ын dL ee AA uer d dapi 75 ph aii dapi өөө эн ад ы дулаан 
2 d e p iiie л АН m A на ne weer 


Два eth inicie EE Aere d a e A Йа У 
ие 
Е Т MEE MR el = 
І > м r (unen A. Pg У 
a d мі га # ыы ӨЛ 2 Та Гоу ја зар Я 
м РР Р С ОО А даран e d рам и 
ЧИ ee diet" к. % 
"пре РУ О pee МР Ања 
i ©» А ө ЭУ ru tet t ova i aa 


^ 





О А 
О AIL) D dee A dai rr rr A b. у ) ^ 
Weird а ми ti e eei dieti АЖ E и : Я . E Eod Ia RA he M па 
Ee eee cha ead hh eee Сы а ој а n" em 8 хал P Ya mck pee аы wh wm ; i à t 
ат аи па ИМ P dad o д М ићи шэн 1 малч лд Л d R T аА ды? ЖЕЛ чом s pum aw » 
bns rdiet Adr de rpm Fede i : ^ к t де Т а њи "DEC. фр а aT Хо, 1 ғ” T қар кт т» 
Г ahn мм МОЛА РА ы; ЕУ У ТА А“ ee “ күк үтү" Ж MEC АКУН 4 LI PLN A ТҮ ТАУДЫ АА ы ЖІ wn was na bell E 
fedi ih f rif APP СР АТА ef, ovS т ES LR d m POPE be AR I Aer OT О ДА ои тоа РЧР РМ И А ы 
2 рт k Ра і = 1 Әм” Тл ТЫ,» МУ ТА АС ТАҚ ы Еу Sut ws БТ Фад ај а Pd њи а АА + 
А А d aw ters 4 ло У АР ЛАРА ПЕК рам МАРАЗ tiep die eL raid qe Lg P NM dq ir ont они niet белье лед ; 
$492 VP Mts cr tla мба ад. 4 ҒЫС”. Талал. dii ут та ето 
^ m Дауни nc bein арад В EE Mid e e арх 
нів a a ӨМ ri eri ra D aaa Aae T 
bbs" O^ St ph арафа і 
тэг” эрй. ти 


